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A long-range interaction via virtual particle-hole pairs between Fermi-liquid quasiparticles leads 
to a nonanalytic behavior of the spin susceptibility x as a function of the temperature (T), magnetic 
field (B), and wavenumber. In this paper, we study the effect of the Rashba spin-orbit interaction 
(SOI) on the nonanalytic behavior of x for a- two-dimensional electron liquid. Although the SOI 
breaks the SU{2) symmetry, it does not eliminate nonanalyticity but rather makes it anisotropic: 
while the linear scaling of x^^ with T and |B| saturates at the energy scale set by the SOI, that 
of Xxx (= Xyy) continues through this energy scale, until renormalization of the electron-electron 
interaction in the Cooper channel becomes important. We show that the Renormalization Group 
fiow in the Cooper channel has a non-trivial fixed point, and study the consequences of this fixed 
point for the nonanalytic behavior of x- An immediate consequence of SOI-induced anisotropy in 
the nonanalytic behavior of x is a possible instability of a second-order ferromagnetic quantum 
phase transition with respect to a first-order transition to an XY ferromagnetic state. 

PACS numbers: 71.10.Ay, 71.10.Pm, 75.40.Cx 



I. INTRODUCTION 

The issue of nonanalytic corrections to the Fermi liq- 
uid theory has been studied extensively in recent yearsjil 
The interest to this subject is stimulated by a vari- 
ety of topics, from intrinsic instabilities of ferromag- 
netic quantum phase transitionsir— to enhancement of 
the indirect exchange interaction between nuclear spins 
in semiconductor heterostructures with potential appli- 
cations in quantum computing Jr^ The origin of the non- 
analytic behavior can be traced to an effective long- 
range interaction of fermions via virtual particle-hole 
pairs with small energies and with momenta which are 
either small (compared to the Fermi momentum kp) or 
near 2kpM-'^ In 2D, this interaction leads to a linear 
scaling of x with a characteristic energy scale E, set 
by either the temperature T or the magnetic field |B|, 
or else by the wavenumber of an inhomogeneous mag- 
netic field |q| (whichever is larger when measured in ap- 
propriate units) J^^— Higher-order scattering processes 
in the Cooper (particle-particle) channel result in addi- 
tional logarithmic renormalization of the result: at the 
lowest energies, x c< E/ In^ E i^^t^^i^^"— The sign of the 
effect, i.e., whether x increases or decreases with E, turns 
out to be non-universal, at least in a generic Fermi liq- 
uid regime, i.e., away from the ferromagnetic instability: 
while the second-order perturbation theory predicts that 
X increases with E, the sign of the effect may be reversed 
either due to a proximity to the Kohn-Luttinger super- 
conducting instability®'^'^^ or higher-order processes in 
the particle-hole channeli^i^ 

In this paper, we explore the effect of the spin-orbit in- 
teraction (SOI) on the nonanalytic behavior of the spin 
susceptibility. The SOI is important for practically all 
systems of current interest at low enough energies; at 



the same time, the nonanalytic behavior is also an inher- 
ently low-energy phenomenon. A natural question to ask 
is: what is the interplay between these two low-energy 
effects? At first sight, the SOI should regularize non- 
analyticities at energy scales below the scale set by this 
coupling. (For a Rashba-type SOI, the relevant scale is 
given by the product |a|fcF, where a is the coupling con- 
stant of the Rashba Hamiltonian; here and in the rest 
of the paper, we set h and fcs to unity). Indeed, as we 
have already mentioned, the origin of the nonanalyticity 
is the long-range effective interaction originating from the 
singularities in the particle-hole polarization bubble. If, 
for instance, the temperature is the largest scale in the 
problem, these singularities are smeared by the tempera- 
ture with an ensuing nonanalytic dependence of x on T. 
On the other hand, if the Zeeman energy /^_b|B| is larger 
than T, it provides a more efficient mechanism of regu- 
larization of the singularities and, as result, x exhibits a 
nonanalytic dependence on |B|. The same argument ap- 
plies also to the |q| dependence. It is often said that the 
SOI plays the role of an effective magnetic field, which 
acts on electron spins. If so, one should expect, for in- 
stance, a duality between the T and \a\kF scalings of x 
(by analogy to a duality between T and /z_b|B| scalings 
of x), be., in a system with fixed SOI, the nonanalytic T 
dependence of x should saturate at T ~ lajfci?. The main 
message of this paper is that such a duality does not, in 
fact, exist; more precisely, not all components of the sus- 
ceptibility tensor exhibit the duality. In particular, the 
in-plane component of Xi Xxx, continues to scale linearly 
with T and yus |B| , even if these energies are smaller than 
lalfciT-. On the other hand, the T and /^bIBJ dependences 
of Xzz do saturate at T ~ |a|/ci^. 

The reason for such a behavior is that although the SOI 
does play a role of the effective magnetic field, this field 
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FIG. 1. (COLOR ONLINE) Rashba spectrum in zero mag- 
netic field. The energy is measured in units of eo = ma^ /2, 
the momentum is measured in units of ma. 



depends on the electron momentum. To understand the 
importance of this fact, we consider the Rashba Hamil- 
tonian in the presence of an external magnetic field ^^ 
which couples only to the spins of the electrons but not 
to their orbital degrees of freedom 
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(1.1) 



where a is the SOI, k is the electron momentum in the 
plane of a two-dimensional electron gas (2DEG), cr is a 
vector of Pauli matrices, e^ is a normal to the plane, g 
is the gyromagnetic ratio, ^b is the Bohr magneton, and 
the effective Rashba field, defined as B^; = (2a/(7/^B)(kx 
e^), is always in the 2DEG plane. The effective Zeeman 
energy is determined by the total magnetic field Btot = 
Bij -h B as 
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(1.2) 



where we introduced A — gnB^/2 for the "Zeeman 
field" , such that A is the Zeeman energy of an electron 
spin in the external magnetic field (h = 1) and 2A equals 
to the Zeeman splitting between spin-up and spin-down 
states. A combined effect of the Rashba and external 
magnetic fields gives rise to two branches of the electron 
spectrum (see Fig. [T]) with dispersions 
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If B II 62, the external and effective magnetic fields 
are perpendicular to each other, as shown in Fig. [2^, 
so that the magnitude of the t otal magnetic field is 
I Btot I = |B -K BflJ = ^52 + B\. This means that the 
external and magnetic fields are totally interchangeable, 
and the T dependence of the spin susceptibility is cut off 
by the largest of the two scales. However, if the external 
field is in the plane (and defines the x axis in Fig. Wp), 




FIG. 2. Interplay between the external magnetic field B and 
effective magnetic field B/j due to the Rashba spin-orbit inter- 
action. Left: the external field is perpendicular to the plane 
of motion. Right: the external field is in the plane of motion. 



the magnitude of the total field depends on the angle 6'k 
between k and B. In particular, for a weak external field. 



Ak ~ |a| A: -I- A sin 6'k. 



(1.4) 



If the electron-electron interaction is weak, the nonana- 
lytic behavior of the spin susceptibility is due to particle- 
hole pairs with total momentum near 2kF, formed by 
electron and holes moving in almost opposite directions. 
In this case, the second term in Eq. (|1.4p is of the op- 
posite sign for electrons and holes. The effective Zee- 
man energy of the whole pair, formed by fermions from 
Rashba branches s and s' is 
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(s-s') |a|fc+ (s + s') A sin 6'k. 

(1.5) 
Only those pairs which "know" about the external mag- 
netic field -via the second term in Eq. (II. 5^ - renormalize 
the spin susceptibility. According to Eq. ()1.5|) , such pairs 
are formed by fermions from the same Rashba branch 
(s = s'). However, since the first term in Eq. p.5p van- 
ishes in this case, such pairs do not "know " about the 
SOI, which means that the Rashba and external mag- 
netic fields are not interchangeable, and the SOI energy 
scale does not provide a cutoff for the T dependence of 

X- 

We now briefiy summarize the main results of the pa- 
per. We limit our consideration to the 2D case and to 
the Rashba SO coupling, present in any 2D system with 
broken symmetry with respect to reversal of the normal 
to the plane. We focus on the dependencies of x o'^ T, 
a, and |i?|, deferring a detailed discussion of the depen- 
dence of X on |q| to another occasion.— Throughout the 
paper, we assume that the SO coupling and electron- 
electron interaction, characterized by the coupling con- 
stant U, are weak, i.e., |a| kp <^ Ep and m[/ ^ 1. The 
latter condition implies that only 2kp scattering pro- 
cesses are relevant for the nonanalytic behavior of the 
spin susceptibility. Accordingly, U is the 2kF Fourier 
transform of the interaction potential. Renormaliza- 
tion provides another small energy scale at which the 
product (r7T,{7/27r) ln(A/rc') (where A is an ultraviolet 
cutoff of the theory) becomes comparable to unity, i.e, 
Tc = Aexp (— 27r/mC/) . Depending on the ratio of the 
two small energy scales-|a| kp and Tc— different behav- 
iors are possible. 
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FIG. 3. (COLOR ONLINE) A sketch of the tempera- 
ture dependence of the spin susceptibihty for the transverse 
(top) and in-plane (bottom) magnetic field. Dashed seg- 
ments in crossovers between various asymptotic regimes do 
not represent results of actual calculations. The left (right) 
panel is valid for Tc <Si \a\kF {Tc ^ |a|fcF) where Tc = 
Aexp (— 27r/mf/) is the temperature below which renormal- 
ization in the Cooper channel becomes significant. The zero 
temperature limit for interacting electrons, denoted by xo, is 
given by Eq. 1)2. 8|) in the Random Phase Approximation. 



In Fig. [3^, we sketch the T dependence of Xzz for the 
case oi Tc <S^ \a\kp. For T ^ \a\kF, Xzz scales linearly 
with r, in agreement with previous studies; a correction 
due to the SOI is on the order of (akp/T) [cf. Eq. 
(|3.29p ]. For Tc 't: T <t: \a\ kp, Xzz saturates at a value 
proportional to |a| ; the correction due to finite T is on 

the order of (T/ |a| kpf [cf. Eq. (l330l) ]. For T < Tc, 
renormalization in the Cooper channel becomes impor- 
tant. In the absence of the SOI, the coupling constant of 
the electron-electron interaction in the Cooper channel 
flows to zero as C//| InTj. Consequently, the T scaling of 
the spin susceptibility changes to T/\v?T for T <^Tc- 
In the presence of the SOI, the situation is different. We 
show that the Renormalization Group (RG) flow of U 
in this case has a non-trivial fixed point characterized 
by finite value of the electron-electron coupling, which 
is only numerically smaller than its bare value. In be- 
tween these two limits, the coupling constant changes 
non-monotonically with InT, and so does Xzz- In both 
high-and low T limits (compared to Tc), however, Xzz is 
almost T independent, so Cooper renormalization affects 
the |a| term in Xzz- 

The T dependence of Xzz for Tc ^ |a| kp is sketched 
in Fig. [3]d. In this case, the crossover between the T and 
T/ \n T forms occurs first, at T ~ Tc, while the T/ In^ T 
form crosses over to |a| at T ^ |a| fcj?. 

We now turn to Xxx- For Tc ^ |a| /cf, its T depen- 
dence is shown in Fig. [3J:. The high-T behavior is again 
linear [also with a (akp/T) correction, cf. Eq. (I3.50p ]. 
For Tc ^ T ^ |a|A:i?, this behavior changes to Xxx 



oc |a| kp/6 + T/2, which means that Xxx continues to de- 
crease with T with a slope half of that at higher T [cf. 
Eq. p.5ip ]. Finally, for T <Tc, Cooper renormalization 
leads to the same |a| dependence as for Xzz but with a 
different prefactor. The behavior of Xxx for |a| kp <^ Tc 
is shown in Fig. [3Ji. Apart from the numbers, this be- 
havior is similar to that of Xzz in that case. 

The dependences of Xxx and Xzz on the external mag- 
netic field can be obtained (up to a numerical coeffi- 
cient) simply by replacing T by A in all formulas pre- 
sented above. In particular, if Cooper renormalization 
can be ignored, both Xxx and Xzz scale linearly with A 
for A 3> |a|A:i? but only Xxx continues to scale with A 
for A <C |a|A:F- A linear scaling of xu with A implies the 
presence of a nonanalytic, \M.i\'^ term in the free energy, 
where M is the magnetization and i — x,y, z. Conse- 
quently, while the cubic term is isotropic {F ex jMl"^) 
for larger M (so that the corresponding Zeeman en- 
ergy is above |Q;|fci?), it is anisotropic at smaller |M|: 
F ex [Mj:!'^ -I- jMj,!'^. A negative cubic term in F im- 
plies metamagnetism and an instability of the second- 
order ferromagnetic quantum phase transition toward 
a first-order one.i-i An anisotropic cubic term implies 
anisotropic metamagnetism, i.e., a phase transition in a 
finite magnetic field, if it is applied along the plane of mo- 
tion, but no transition for a perpendicular field, and also 
that the first-order transition is into an XY rather than 
Heisenberg ferromagnetic state. This issue is discussed 
in more detail in Sec. |Vl 

The rest of the paper is organized as follows. In SecHIl 
we formulate the problem and discuss the T dependence 
of the spin susceptibility for free electrons with Rashba 
spectrum (more details on this subject are given in Ap- 
pendix |X|. Section IlII Al explains the general strategy of 
extracting the nonanalytic behavior of x from the ther- 
modynamic potential in the presence of the SOI. The 
second-order perturbation theory for the temperature 
and magnetic-field dependences of the transverse and in- 
plane susceptibilities is presented in Sees. IIIIBI and llllCl 
respectively. In Sec. IIII Dl we show that, as is also the 
case in the absence of the SOI, there is no contribution to 
the nonanalytic behavior of the spin susceptibility from 
processes with small momentum transfers, including the 
transfers commensurate with (small) Rashba splitting of 
the free-electron spectrum (more details on this issue are 
provided in Appendix IB)) . Renormalization of spin sus- 
ceptibility in the Cooper channel is considered in Sec.lIVI 
An explicit calculation of the third-order Cooper contri- 
bution to Xzz is shown in Sec. IIVBI In Sec. lIVC^ we de- 
rive the RG fiow equations for the scattering amplitudes 
in the absence of the magnetic field; the effect of the fi- 
nite field on the RG fiow is discussed in Appendix [C] The 
effect of Cooper-channel renormalization on the nonana- 
lytic behavior of Xzz and Xxx is discussed in Secs. lIVCSI 
and IIV C 41 correspondingly. Implications of our results 
in the context of quantum phase transitions are discussed 
in Ref. |Vl where we also give our conclusions. 



II. SPIN SUSCEPTIBILITY OF FREE RASHBA 
FERMIONS 

In this Section, we set the notations and discuss briefly 
the properties of Rashba electrons in the absence of the 
electron-electron interaction. The Hamiltonian describ- 
ing a two-dimensional electron gas (2DEG) in the pres- 
ence of a Rashba SOI and an external magnetic field B 
is given by Eq. (|l.ip . 

In the following, we consider two orientations of the 
magnetic field: transverse (B = Be^) and parallel (B — 
Be^) to the 2DEG plane. It is important to emphasize 
that, when discussing the perpendicular magnetic field, 
we neglect its orbital effect. Certainly, if the spin suscep- 
tibility is measured as a response to an external magnetic 
field, its orbital and spin effects cannot be separated. 
However, there are situations when the spin part of Xzz 
is of primary importance. For example, the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interaction between the 
local moments located in the 2DEG plane arises only 
from the spin susceptibility of itinerant electrons, be- 
cause the orbital effect of the dipolar magnetic field of 
such moments is negligible. In this case, Xxx and Xzz de- 
termine the strength of the RKKY interaction between 
two moments aligned along the x and z axis, respec- 
tively. Also, divergences of Xzz and Xxx signal ferromag- 
netic transitions into states with easy-axis and easy-plane 
anisotropics, respectively. Since it is this kind of physical 
situations we are primarily interested in this paper, we 
will ignore the orbital effect of the field from now on. 

The Green's function corresponding to the Hamilto- 
nian (jl.ip is obtained by matrix inversion 
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where we use the "relativistic" notation K = (w, k) with 
cj being a fermionic Matsubara frequency, the matrix 
ris(k) is defined as 

1 
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0«(k) ; 
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a{kyd-x 



l^x^y) 
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(2.2a) 
(2.2b) 



gs{K) — l/(ia; — ^k — sAk) is the single-electron Green's 

function, ^k = j /x, /i is the chemical potential, and 

Ak is given by Eq. (fr2l) . 

As we have already pointed out in the Introduction, an 
important difference between the cases of transverse and 
in-plane magnetic field is the dependence of the effective 
Zeeman energy [Eq. (|1.2p ] on the electron momentum. 
For the transverse magnetic field, (A^; = Aj^ = 0, Az = 
A), the Zeeman energy is isotropic in the momentum 
space and quadratic in A in the weak-field limit: Ak ~ 
\a\kp 4- A^/2|a|fci?. Correspondingly, the Fermi surfaces 
of Rashba branches are concentric circles with slightly 
(in proportion to A^) different radii. For the in-plane 
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FIG. 4. (COLOR ONLINE) Spin susceptibility of free Rashba 
fermions (in units of xo = Mb'ti/tt) as a function of the 
chemical potential /i (in units of eo — ma^ /2). Solid (red) 
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X^S ; dashed (blue): x 



Note that xi? 



/ Xzz if only 



the lowest Rashba branch is occupied (— eo < M < 0) but 
Xxx — Xzz — Xo if both branches are occupied {fi > 0) at 

r = o. 



Zeeman energy is anisotropic in the momentum space, 
cf. Eq. (|1.4p . Correspondingly, the Fermi surfaces of 
Rashba branches are also anisotropic and their centers 
are shifted by finite momentum, proportional to A^;. 

We now give a brief summary of results for the sus- 
ceptibility in the absence of electron-electron interaction, 
which sets the zeroth order of the perturbation theory 
(for more details, see Appendix \Kl . The in-plane rota- 
tional symmetry of the Rashba Hamiltonian guarantees 
that Xvv — Xxx ■ The static uniform susceptibility (de- 
fined in the limit of zero frequency and vanishingly small 
wavenumber) is still diagonal even in the presence of the 
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general. The susceptibility depends strongly on whether 
both or only the lower of the two Rashba branches are oc- 
cupied, see Fig. (j4]). In the latter case, the spin response 
is strongly anisotropic. At T = 0, 



xf) == XoVl + M/eo 
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A,XX 
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(2.3a) 
(2.3b) 



where xa = [^b^^I'^ i^ ^^^ spin susceptibility of 2D elec- 
trons in the absence of the SOI, eo = mo? jl is the depth 
of the energy minimum of the lower branch and the chem- 
ical potential [i is within the range — Eq < /i < 0. The in- 
plane susceptibility exhibits a ID-like van Hove singular- 
ity at the bottom of the lower branch, i.e., for /i — > — eo. 
On the other hand, if both branches are occupied (which 
is the case for /i > 0), the spin susceptibility is isotropic 
and the same as in the absence of the SOI 
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(2.4) 



magnetic field, (A^; — A, Aj, = A^ = 0), the effective This isotropy can be related to a hidden symmetry of 



+ <a:;> + <ri> + 



FIG. 5. The RPA diagrams for the spin susceptibility corre- 
sponding to Eq. fTf}. 



the Rashba Hamiltonian manifested by conservation of 
the square of the electron's velocity operator v^ The 
eigenvalue of w^ , given by 2e/m + 2a^ with e being the 
energy, is the same for both branches. The square of the 
group velocity v^ — (Vkf^ ) — 2e/m + a^ also does not 
depend on the branch index. Therefore, the total density 
of states 
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where k = ^ma + \/vn?-c? + 2me are the momenta of 
I 



the ± branches corresponding to energy e, is the same 
as without the SOI, if both branches are occupied. One 
can show also that isotropy of the spin susceptibility is 
not specific for the Rashba coupling but is there also 
in the presence of both Rashba and (linear) Dresselhaus 
interactions;^ 

As one step beyond the free-electron model, we con- 
sider a Stoner-like enhancement of the spin susceptibility 
by the electron-electron interaction. In the absence of the 
SOI and for a point-like interaction [/, the renormalized 
spin susceptibility is given by 
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In the presence of the SOI, the ladder series for the sus- 
ceptibility, shown in Fig. [5] is given by 
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where Q = (fi = 0,q — > 0) and i — x^z. Although the 
traces do look different for Xxx and Xzz , these differences 
disappear after angular integrations, and the resulting 
series are the same. As we have shown above, the zero- 
order susceptibilities are also the same (and equal to Xo) 
if both Rashba branches are occupied; hence 
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(2.8) 
which is the same result as in Eq. (|2.6p . Therefore, at the 
mean-field level, the spin susceptibility remains isotropic 
and independent of the SOI. In the rest of the paper, 
we will show that none of these two features survives 
beyond the mean-field level: the actual spin susceptibility 
is anisotropic and both its components do depend on the 
SOI. 

We now come back to the free-electron model and dis- 
cuss the T dependence of the spin susceptibility. A spe- 
cial feature of a 2DEG in the absence of the SOI is a 
breakdown of the Sommerfeld expansion at finite T: since 
the density of states docs not depend on the energy, all 
power-law terms of this expansion vanish, and the result- 
ing T dependence of x^^^ is only exponential. The SOI 
leads to the energy dependence of the density of states 
for the individual branches, and one would expect the 
Sommerfeld expansion to be restored. This is what in- 
deed happens if only the lower branch is occupied. In 
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{0}, 



provided that T <C min {— ^, eq + A^}- [Here, xL (0) and 
xix (0) are the zero temperature values given by Eqs. 
(I2.3a|) and (I2.3bp ]. However, if both branches are occu- 
pied, the energy dependent terms in the branch densities 
of states cancel out, and the resulting dependence is ex- 
ponential, similar to the case of no SOI, although with 
different pre-exponential factors: 



x|°>(T)=xo(l-U^e--/^ 
xS>(T)=Xo(l + ^e-''/^ 



for T ^ £0 ^ /^, and 
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(2.10a) 
(2.10b) 

(2.11a) 
(2.11b) 



for eo <C T <C /x. In a similar way, one can show that 
that there are no power-law terms in the dependence of 

Xxx and Xzz on the external magnetic field. 

Notice that Xxx 7^ Xzz at finite temperature, even if 
both Rashba subbands are occupied. This suggests that 
the hidden symmetry of the Rashba Hamiltonian is, in 
fact, a rotational symmetry in a (2 + 1) space with imag- 
inary time being an extra dimension. ^^ Finite tempera- 
ture should then play a role of finite size along the time 
axis breaking the rotational symmetry. 

In what follows, we assume that the SOI is weak, i.e., 
I a I <^ vf, and thus the energy scales describing SOI are 



small: 



<^ \a\kF <^ 11. This condition also means 



that both Rashba branches are occupied. [Also, from now 
on we relabel fj, — >■ Ep-] The main result of this section is 
that, for a weak SO coupling, the T and A dependences of 
X in the free case are at least exponentially weak and thus 
cannot mask the power-law dependences arising from the 
electron-electron interaction, which are discussed in the 
rest of the paper. 



To second order in the electron-electron interaction U{q), 
there is only one diagram for the thermodynamic poten- 
tial that gives rise to a nonanalytic behavior: diagram 
a) in Fig. [51 The rest of the diagrams in this figure can 
be shown to be irrelevant (cf. Sec lIIIDJ) . Algebraically, 
diagram a) in Fig. [6] reads 



^-^'^--iEEE'^'(ik-pi) 



Q K P 

X Tr{GKGp)Tr{GK+QGp+Q), (3.2) 

where ^q = Tj:ald'q/{27r)^ Ek = 
T^j^ /d^fc/(27r)^, and we use "relativistic" nota- 
tion K = (w, k) with a fcrmionic frequency uj and 
Q = (ri, q) with a bosonic frequency f2. 

Evaluation of the first spin trace in Eq. p.2p yields 



III. SPIN SUSCEPTIBILITY TO SECOND 

ORDER IN THE ELECTRON-ELECTRON 

INTERACTION 
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A. General strategy 



where 



The spin susceptibility tensor Xij is related to the ther- 
modynamic (grand canonical) potential S(T, a, A) by the Bss' (k, p) = 1 -|- ss' 
following identity 



, a^k ■ p + g [A X (k + p)]^ + A^ 
AkAp 



Xzj (T, a) 



d^E 



dBidBj 



B=0 



(3.4) 

and Ak is given by Eq. (|1.2p . The second-order thermo- 
(3.1) dynamic potential then becomes 

I 



^S(2) = -^ EE E ^'(\^ - p\)Bs,s,{Kp)Bs,sA^ + q,P + ci)gsdK)9s,{P)9sdK + Q)9sAP + Q)- (3-5) 



K P 



Equation p.5|) describes the interaction among electrons 
from all Rashba branches via an effective vertex [/(|k — 
p\)Bss' (k, p), which depends not only on the momentum 
transfer k — p but also on the initial momenta k and p 
themselves. This last dependence is due to anisotropy of 
the Rashba spinors. 

It has been shown in Refs.y and [la that, at weak cou- 
pling, the main contribution to the nonanalytic part of 
the spin susceptibility comes from "backscattering" pro- 
cesses, i.e., processes with p w — k and small q (compared 
to kp)- In particular, the second-order contribution is 
entirely of the backscattering type. The proof given in 
Ref. [3] applies to any kind of the angular-dependent ver- 
tex and, thus, also to vertices in Eq. p.Sp . Therefore, 
the calculation can be simplified dramatically by putting 
in p = — k and neglecting q in the effective vertices. The 
last assumption is justified as long as the typical values of 



r 



q are determined by the smallest energy of the problem, 
i.e., q ^ max{T/wi?,r7T,|a|} <C kp- By the same argu- 
ment, the magnitudes of k and p in the vertices can be 
replaced by kp. The bare interaction is then evaluated 
at |k — p| « 2kp, and we introduce a coupling constant 
U = U{2kp). Ignoring the angular dependence of Bs^s 
for a moment, Eq. (j3.5p is reduced to a convolution of 
two particle-hole bubbles, formed by electrons belonging 
to either the same or different Rashba branches 



ns..,(0) = 5]<?.,(i^)ff.,(if + Q). 



(3.6) 



K 



By assumption, both components of Q in Ils^sjiQ), i-©-, 
fl and g, are small (compared to Ep and kp, correspond- 
ingly). It is important to realize that, despite a two- 
band nature of the Rashba spectrum, Hs^s (Q) has no 
threshold-like singularities at the momentum qo = 2ma, 



U{k-p 




U{k - p) [;(,/) 




i) 



U{q) K - Q 




K + Q 




K + Q' 



K + Q 




K + Q' 



FIG. 6. Second order diagrams for the thermodynamic potential. A nonanalytic contribution comes only from diagram a), 
where the momenta are arranged in such a way that k ~ — p while Q is small; therefore, k ^ p ^ kp, and the momentum 
transfer in each scattering event is close to |k — p| « 2kF- As is shown in Sec. IIIID] diagrams 6)-e) do not contribute to the 
nonanalytic behavior of the spin susceptibility. 



separating the Rashba subbands2^ (for a detailed deriva- 
tion of this result, see Appendix IB]) . Therefore, the non- 
analytic behavior of the spin susceptibility comes only 
from the Landau-damping singularity of the dynamic 



bubble, as it is also the case in the absence of the SOI. 



After the simplifications described above, Eq. 
comes 



3be- 



U 



SE^''> = -—Y,J2J2^^^^^^^^-^-9sAK)gsAP)9sdK + Q)9sAP + Q), 



(3.7) 



where 

Bss' = Bss' i'kp,~'kp) = 1 -f- ss' 



A2 



a'^kp 



Ak.A_i 



(3.8) 



and ki^ = {\i./k)kp. Finally, to obtain the leading T 
dependence of Xiji it suffices to replace the fermionic 
Matsubara sums in Eq. p.7|) by integrals but keep the 
bosonic Matsubara sum as it is. The rest of the calcu- 
lations proceed somewhat differently for the cases of the 
transverse and in-plane magnetic fields. 



B. Transverse magnetic field 

First, we consider a simpler case of the magnetic field 
transverse to the 2DEG plane: A = gusBez/'i,. In this 
case, the effective Zeeman energy is isotropic in the mo- 
mentum space; therefore, Ak^ = A_kjj and 



B^f,' = 1 



,A^-a^fc| 



1 + ss' 



A^ + a'^kj,' 



(3.9) 



Thereby the integrals over d^K and cfiP separate, and 
one obtains 



(5S(2) 



El 
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/ ^ -^SlS3-^S2S4 / ^ ^^SiS2W)^^S3SaW)^ 



S1...S4 — ±1 



(3.10) 

where IIs-s (Q) is given by Eq. p.6p . The single-particle 
spectrum in the second Green's function in Eq. p.6p can 
be linearized with respect to q as ^k+q ~ ^k + VFQ cos (pkq 
with (/)kq = Z(k, q). Since, by assumption, the SOI is 



I 

small, the effective Zeeman energy in the Green's func- 
tions can be replaced by its value at k — kp. Integration 
over dkk can be then replaced by that over d^k- The in- 
tegrals over w, ^k and (pkq are performed in the same way 
as in the absence of the SOI, and we arrive at the follow- 
ing expression for the dynamic part of the polarization 
bubble 



n. 



101 



2tt 



"n + ^(s'~s)Ak^] +{vpqf 



(3.11) 



We pause here for a comment. The polarization bubble 
in Eq. p.ll|) is very similar to the dynamic part of the 
polarization bubble for spin-up and -down electrons in 
the presence of the magnetic field but in the absence of 
the SOI^ 



nn(o,<z) = 



101 



(3.12) 



2n ^{n - igfiBB)^ + [vpqf 

As we already mentioned in Sec. HI the nonanalytic be- 
havior of the spin susceptibility is due to an effective in- 
teraction of Fermi-liquid quasiparticles via particle-hole 
pairs with small energies and momenta near 2kp. Our 
calculation is arranged in such a way that, on a tech- 
nical level, we deal with pairs with small momenta q. 
The spectral weight of these pairs is proportional to the 
polarization bubble, which is singular for small O and 
q. Since finite magnetic field cuts off the singularity in 
n-|-4,(0, g), the nonanalytic dependence of x on, e.g., tem- 
perature, saturates when T becomes comparable to the 
Zeeman splitting g^sB. At lower energies, x exhibits 
a nonanalytic dependence on A: x oc |A|. Likewise, the 
singularity in Eq. p. lip is cut at the effective Zeeman en- 
ergy Akp, which reduces to the SOI energy scale lajfej^. 



when the real magnetic field goes to zero. Therefore, one 
should expect the nonanalytic T dependence of Xzz to 
be cut by the SOI. Although soft particle-hole pairs can 
be still generated within a given branch, i.e., for s = s' ^ 
the entire dependence on the Zeeman energy in this case 
is eliminated and processes of this type do not affect the 
spin susceptibility. In the rest of this section, we are go- 
ing to demonstrate that the SOI indeed plays a role of 
the magnetic field for Xzz- 

For later convenience, we define a new quantity 



Vs^Vs{n,q) 
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[n - 2isAkj 



(3.13) 
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Vpq 



and sum over the Rashba branches in Eq. (|3.10l) . The 
contribution of the set {si = S2 = s, S3 = S4 = s'} does 
not depend on A\^p : 



n 



E^i'nLH^ E^i'^o^ 



(3.14) 



and obtain, after some algebra, 

{STrJ ^ Jo 27rA4/ ^- 

-FA''(p2 +vl- 2Vi) + Aa^kj, A:^Vo{V+ +V-- 2750)], 

(3.18) 
where we used that / dqq {V+V- - Vl) = 0. The first 
term in the square brackets in Eq. p.lSp does not depend 
on the effective field Akj^ and, therefore, can be dropped. 
Integration over q in the remaining terms is performed as 



/ 



dqqVo{V+ +V-- 2Vo) = ^ In ^ " . , , (3.19) 



■jI n^ + Kl. 



( dqq{Vl+Vl-2Vl) = \\n 
J Vp 



n^ 



n^ + 4A2 



(3.20) 



Collecting all terms together, we obtain 



The set {si = —S2 = s, S3 = — S4 = s'} gives 



/ ^ Bss' B-s,-s'^s-s^s 



rnfl 



V HBl'V.V. 



s,s' 



(3.15) 
where we used that B^g = B_s,-s- Finally, the sets {si = 
S2 = s, S3 = -S4 = s'} and {si = -S2 = s, S3 = S4 = s'} 
contribute 



TT \ AtTVF J 
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(3.21) 



The bosonic sum is evaluated by replacing the sum by an 
integral as follows 



y{Bss'Bs,^s,UgsUs,^^s'+Bss,B^s,s>n,^^sns's') TYF(n)^ H — coth— Im 



lim F(~in + S) 



= 2(^yj2Bgg,Bg^_g,roVs,, 

^ '^ ss' 

(3.16) 

where we relabeled the indices in the second sum {s -^ s' ^ 
s' — > s) and used the symmetry property Egsi — Ba's- 
The angular integral contributes a unity, so that 






fjiU 



SW ^ 7 2vr 

ss' 

(3.17) 



Now we sum over ss'^ add and subtract a combination 
2(Q:^fc|, -I- ia'^kpA'^ + A'')7'o inside the square brackets. 



(3.22) 



and using the identity 

Im [ lim {-2n + 5f In^i-^^^-tiL 

= TTVL'^s\gnnQ{x^ -9^), (3.23) 

where Q{x) stands for the step function. 

The thermodynamic potential then becomes 
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(3.24) 



with 



T{y) = / dxx' coth (.t/2) = -{l/3)y' [y - 
Jo 



6Lii (e^)] + 4[C(3) + yU2 (e^) - Lis (e^)], 



(3.25) 



where Li„(z) = X]fe=i z /k" is the polylogarithm func- tion and C(-2) is the Riemann zeta function. In practice. 



the integral form of T{y) is more convenient as it can 
be easily expanded in the limits of small and large ar- 
gument. Indeed, for j/ ^ 1, one expands coth(a;/2) as 
coth(a;/2) = 2/x + x/6 and, upon integrating over x in 



Eq. (|3^ . obtains 



•^(y) = y' + f^ + o (2/') , for < y « 1. 



(3.26) 



For y 3> 1, one subtracts unity from the integrand and 
replaces the upper limit in the remaining integral by in- 
finity: 



:f{v) 



3 



dxx I coth 



1 W ^ + / dxx"^ f coth I - 1) - / dxx'^ f coth ^ - 1 



4C(3) + 0(e' 



for jy >• 1. 



(3.27) 



1. Temperature dependence of Xzz 

The (linear) spin susceptibility is given by Xzz = 
—d'^'E,/dB^\B=o, which means that only terms propor- 
tional to A^ in the thermodynamic potential matter. 
Therefore, for finite a, the spin susceptibility comes en- 
tirely from the second term in the square brackets of 
Eq. (|3.24l) . which is proportional to a^fc|,A^. [On the 
other hand, for a — the second term vanishes, while 
A'^/A^ = 1, and the spin susceptibility comes exclu- 

sively from the first term: S^zz still depends on the mag- 
netic field through J^(Akjj/T), where the A-dependence 
must be retained; in this case, J^(Akj^/T) is evaluated 
as shown in Ref. [l^.] Neglecting the first term and dif- 
ferentiating the second one, we obtain the interaction 
correction to Xzz for A — > as 



^^^-»(f)^^ 



T 



\a\k[ 
T 



(3.28) 



For T 3> |a|fcF, the asymptotic expansion of J- in Eq. 
(1?:^ gives 



5xf) « 2x0 



mU 
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1 a^kj, 
2ATEf 



(3.29) 



The first term in Eg. (| 3.29l) coincides with the result of 
Refs. y, [1^ [21I, and |22 obtained in the absence of the 
SOI, while the second term is a correction due to the 
finite SOI. In the opposite limit, i.e., for T <C |a|fcF, the 
asymptotic expansion of J- in Eq. (I3.27|) gives 



5x 
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a'^kj.Ep 
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(3.30) 
As it was anticipated, the SOI cuts off the nonanalytic T 
dependence for T < jajfei?- However, the T dependence 
is replaced by a nonanalytic \a\ dependence on the SO 
coupling. 

Normalizing Eq. (|3.28l) to the leading T depen- 
dent term for a = 0, i.e., by 5xzz{T,a — {)) — 
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SSI.5 







0.5 
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T/|a|^/ 



FIG. 7. (COLOR ONLINE) The second-order nonanalytic 
correction to x^^y normalized to its value in the absence of the 
SOI, as a function oiT/\a\kF [cf. Eq. (I3.31[l ]. The asymptotic 
forms, given by Eqs. (|3.29(l and (|3.29p . are shown by dashed 
(red and blue) lines. 



Xo (mC//47r) [T/Ep) , we express 
function of the variable T / \a\kp 



X (2) 

6Xzz 



via a scaling 



5xi^}{T,a) 
SxfhT.a 



0) 



T 



\kf 



T 



\oi\ki 
T 



(3.31) 



The left-hand side of Eq. (|33T|) is plotted in Fig. [7] along 
with its high and low T asymptotic forms. 

Now we can give a physical interpretation of the above 
results. Although the electron-electron interaction mixes 
Rashba branches, the final result in Eq. p.28p comes only 
from a special combination of electron states. Namely, 
three out of four electron states involved in the scattering 
process (two for the incoming and two for the outgoing 
electrons) must belong to the same Rashba branch, while 
the last one must belong to the opposite branch, as shown 
in Fig. [B^. This can be seen from Eq. (|3.I8p by consider- 
ing four terms in the square brackets. The first term, pro- 
portional to A^ , does not depend on the field upon the 
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+ 



b) 





+ 



FIG. 8. Top: Diagrams contributing to nonanalytic behavior 
of x^z- There are eight such diagrams with the following 

choice of Rashba indices: + + + — , + H h, H \ — h, — h H — h, 

h, h — , \ , and H ; one of the them 

is shown on the right. Bottom: Diagrams contributing to 
nonanalytic behavior of Xxx- There are four such diagrams 

with the following choice of Rashba indices: H h — , H h, 

— h + — , and — I h; one of them is shown on the right. 



cancellation with an overall factor of At in the denom- 
inator; the second term, proportional to a'^kp, vanishes; 
the third term is already proportional to A** and, thus, 
cannot affect the spin susceptibility for finite a. There- 
fore, the effect comes exclusively from the last term, pro- 
portional to A^, because the product Bss'Bg^^s' is equal 
to AA^iakpf /{a^kj, + A^)^ « AA^/iakp)^ for any 
choice of s and s' . This term corresponds to the structure 
Y.^^,Bss'Bs,-s'n,^^sUs's' in Eq. (IXT71) . The diagrams 
corresponding to this structure are shown in Fig. |5^. 
Pairing electron Green's functions from different bubbles, 
we always obtain a combination ^^ g± [K) g± {K + Q) , 
which depends neither on the SOI nor on the mag- 
netic field, and a combination J2k 5± i-^) 9t (^ + Q) ' 
which depends on both via the effective Zeeman energy 
Akp = \J a^k\ -f A^. In the weak-field limit, one needs 
to keep A^ only in the prefactor. The singularity in the 
combination "YIiK •9± (^) 5=F (^ ^ Q) ^^ then regularized 
by finite SOI, which is the reason why the nonanalytic T 
dependence is cut off by the SOL 



2. Magnetic-field dependence of Xzz 

Now we consider the case of T <C max{A, |a|fcF} when, 
to first approximation, one can set T = 0. In this case, 
one can define a non-linear susceptibility Xzz iBz,a) — 
—d'^Ezz/dBl evaluated at finite rather than zero mag- 
netic field. For T ^- 0, we replace the scaling function 
J^ in Eq. (j3.24p by the first term in its large- argument 
asymptotic form p.27|) to obtain 
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(3.32) 



Differentiating twice with respect to the field, we find 
«<?<B.») = »(^)^^<;(^), (3.33, 



where 



Six) 



2x^ + 23x^ + 30a;-^ 
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(3.34) 



3(l + a;2)^/2 

has the following asymptotics 

g (x < 1) = 4 + xVe + . . . 

g {x :$> 1) ^ {2/3)\x\ + 6/\x\ + . . . 

For |A| ^ |a| kp, the nonanalytic correction Xzz {Bz,Oi) 
reduces to the result of Ref. |3|, obtained in the absence 
of the SOI, plus a correction term 



(3.35) 



<5xE^(S„a) = 4xo 
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1 a^kj, 



Ef ^aiai-'Ef 



(3.36) 

In the opposite limit of |A| ^ |a| fef, the nonanalytic 
field dependence is cut off by the SOI 
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(3.37) 



C. In-plane magnetic field 

If the magnetic field is along the a;-axis, A = 
fiBBex/2^ the effective Zeeman energies, Aj-k^ = 
y/a'^kp ± 2akF sin6'kA + A^, depends on the angle 0k 
between k and the direction of the field, chosen as the x 
axis. Coming back to Eq. p.7p . we integrate first over 
the fermionic frequencies, then over the magnitudes of 
the fermionic momenta, then over the angle between p 
and q, and finally over the angle between q and k (at 
fixed k). This yields 
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(3.39) 



with gf'^'^ (K) — l/(ia; — ^k — sAikp), and J2k indicates 
that the integration over 9k is excluded. The remaining 
integration over 9k is performed last, after integration 
over q and summation over fl. 
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As in Sec. IIII Bl it is convenient to define a new quantity 



'f}-2JsA±k^) +vW 



and to re-write the thermodynamic potential as 



5^ = - (^) '^E^VIr/lr^ [Bs.B^.^.Poiv:;^^ + v;,^^ ) + Bl, {vl + vr^v; 



(3.40) 



) • (3.41) 



Subsequently, we sum over s and s' , add and subtract 4(A^ A^j, — 2a^fc|, sin20k)'Po inside the square brackets, 
and, after some algebraic manipulations, obtain 



a+ ( T'+^'^T'I*'^ + 7?7kF7?+kF _ 2-p^ ) + a_ ( -p-kF-p+k^ ^ p-k^p+k^ _ ^n 



where 
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4a2fc|,A2 cos^ 6*1^ 



-; a± 



a^k% + A^ - 2a'kj,A' sin^ 6>k ± Ak^A„k^(a^fcj. - A^ 



(3.42) 



(3.43) 



The first term in the square brackets in Eq. (13.421) does not depend on the effective field A±kir and can be dropped. 
The remaining integrals over q are equal to 



JdqqVo[P-^ 
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Evaluating the Matsubara sum in the same way as in Sec. IIIIBi we obtain 
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(3.46) 



(3.47) 



where the function J^ (x) and its asymptotic limits are 
given by Eqs. (I3.25l f 



The angular integral cannot be performed analytically 
because the function J^{y) depends in a complicated way 
on the angle 0k through A±kF- Therefore, we consider 
two limiting cases below. 



1. Temperature dependence of Xxx 

First, we consider the limit of a weak magnetic field: 
|A| <C max{|a|fci?,T}. The main difference between the 
cases of in- and transverse orientations of the field is in 
the term proportional to a+ in Eq. p.47p . The argument 
T in this term vanishes in the limit of A — ^ 0, whereas 
the arguments of T in the rest of the terms reduce to a 
scaling variable \a\kp/T, as it was also the case for the 
transverse field. Therefore, the SOI energy scale, lalfcj-. 
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and temperature are interchangeable in the rest of the 
terms, which means that a nonanalytic T dependence 
arising from these terms is cut off by the SOI (and vice 
versa). However, the 04. term does not depend on a and 
produces a nonanalytic T dependence which is not cut 
off by the SOI. To see this, we expand prefactors ao and 
a± to leading order in A as 
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(3.48) 



The last term in Eq. p.47p . proportional to a_, does 
not contribute to order A'^ , and we focus on the first two 
terms. In the term ao, we replace A±kj? = |a| kp in the 
argument of the J^ function. In the a+ term, we replace 
|Ak^ - A_kf I = 2|Asin6'k| + ^(A^), and th en ex pand 
J"(2|Asin6'k|/r) = iA'^ sin^ Ou/T'^ using Eq. (IXMl) . In- 
tegrating over 9k and differentiating the result twice with 
respect to B, we obtain 
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for T> \a\kF and 
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for T < |a|fcF. 



pEp 
(3.51 



As it is also the case for the transverse magnetic field, 
the first term in the second line of Eq. p.50|) is due to 
particle-hole pairs formed by electrons from three iden- 
tical and one different Rashba branches. On the other 

f2) 

hand, the linear- in-T term, absent in Sxzz , is comes from 
processes involving electrons from the different Rashba 
branches in each particle-hole bubble, see Fig. [5}d. In- 
deed, pairing electrons and holes, which belong to dif- 
ferent Rashba branches and move in the same direction, 
we obtain particle- hole bubbles V^^" [cf. Eq. dHU])]. 
The product of two such bubbles, V^ *"p^ ^^ being in- 
tegrated over q and summed over 17, depends on the dif- 
ference of the Zeeman energies |Akjj— A^kpl- Since sin 0k 
is odd upon k — ?> — k, this difference is finite and propor- 
tional to |A| for A — > but does not depend on a. This 
is a mechanism by which one gets an ©(A^) contribution 
to the thermodynamic potential and, therefore, a T de- 
pendent contribution to Xxx, which does not involve the 
SOI. 



where Xzz {T, a) is the correction to Xzz given by Eq. 
(|3.28|) . [Notice a remarkable similarity between Eq. 
p.49|) and the relation between Xxx and Xzz iu the non- 
interacting case, Eq. (IA5|) .] Equation (|3.49p is one of the 
main results of this paper. It shows that a nonanalytic 
T dependence of Xxx, given by the stand-alone TjEp 
term survives in the presence of the SOI. Explicitly, the 
T dependence is 
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2. Magnetic-field dependence of Xx 



Now we analyze the non-linear in-plane susceptibil- 
ity Xxx{B x,a) = -d^E^a:/dBl at T = 0. Replacing 
J^ in Eq. (j3.47p by its large- argument asymptotics from 
Eq. ((X?fl) . we obtain 



5S(2) 
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8nvp 



i2ny 



{ao[ 



AL 



+ At 



a+|Akj. 



A_ 



+ a- [At 



A- 



'}■ 



(3.52) 



r 



The angular integral can now be solved explicitly in the 
limiting cases of |A| ^ \a\kp and |A| ^ |a|A:_F. Since 
our primary interest is just to see whether a nonanalytic 
field-dependence survives in the presence of the SOI, we 
will consider only the weak-field case: |A| ^ |a|fci?. The 
a_ term in Eq. p.52p can then be dropped, while the oq 



and a+ ones yield 



Sx^xx{Bx,a) = -xo 



mU 

An 



2 r 



\a\ kf 

~b7 



16^ 
n Ep 



(3.53) 
The first term in Eq. p.53|) is just half of the first term 

in Sxzz [cf. Eq. p.37|) ]. However, the second term repre- 
sents a nonanalytic dependence on the field which is not 
cut off by the SOI. 
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D. Remaining second order diagrams 

Besides the diagrams considered so far, there are other 
second order diagrams, which ~ in principle - could con- 
tribute to the spin susceptibility. These diagrams are 
depicted in Fig. [B^-e. In the absence of the SOI, these 
diagrams are irrelevant because the electron-electron in- 
teraction conserves spin. This means that the spins of 
electrons in each of the bubbles in, e.g, diagram b) are 
the same and, therefore, the Zeeman energies, entering 
the Green's functions, can be absorbed into the chemical 
potential. The same argument also goes for the other 
two diagrams. In the presence of the SOI, this argu- 
ment does not work because spin is not a good quantum 
number and the interaction mixes states from all Rashba 
branches with different Zeeman energies. However, one 
can show that the net result is the same as without the 
SOI: diagrams in Fig. |Bf)-e do not contribute to the non- 
analytic behavior of the spin susceptibility. This is what 
we are going to show in this section. 

We begin with diagram 6), which is a small 
momentum-transfer counterpart of diagram in a): 



.Sf) 



-U 



'■{0)TY,[T^rn{Q) 



(3.54) 



where 



n(Q) = ^E^(^)^(^' + Q) 

K 



s,t=±l 



K 



is the full (summed over Rashba branches) polarization 
bubble, and both the space- and time-like components 
of Q = (f^,q) are small. As it is also the case in the 
absence of the SOI, the small-Q bubble does not depend 
on the magnetic field. Indeed, noticing that the matrix 

I 



( in Eq. (|2.2ap has the following properties 

C^ = / and TrC = 0, 
it is easy to show that 
1 



".s"t — — 



{l + st)I + {s + t)C 



(3.56) 



(3.57) 



Consequently, fi+f2+ = 51_ri_ = (1/2)/ and fj+fi- = 
Ct^Cl^ = 0. Therefore, electrons from different branches 
do not contribute to 11 (Q), while the Zeeman energies in 
the contributions from the same branch can be absorbed 
into the chemical potential. As result, the dynamic part 
of the bubble depends neither on the field nor on the SOI 
(as long as a weak dependence of the Fermi velocity for 
a given branch on a is neglected): 



n(Q) 



\n\ 



v^q^ 



(3.58) 



2t: ,/W 
Therefore diagram b) does not contribute to the spin sus- 
ceptibility. We remind the reader that, since there are no 
threshold-like singularities in the static polarization bub- 
ble (see discussion in Sec. IIII| . the Landau-damping sin- 
gularity in Eq. p. 581) is the only singularity which may 
have contributed to a nonanalytic behavior of the spin 
susceptibility. However, as we have just demonstrated. 
Landau damping is not effective in diagrams with small 
momentum transfers. 

Similarly, diagram c) with two crossed interaction 
lines, one of which carries a small momentum and the 
other one carries a momentum near 2kp, is expressed via 
a small Q bubble as 



(5S(2) = _c/(0)C/(2fc, 



.)5]Tr[n2(Q) 
Q 



(3.59) 



and, therefore, does not depend on the magnetic field. 
Diagram d), with both interaction lines carrying small 
momenta, contains a trace of four Green's functions 
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(3.60) 
(3.61) 



I 

where Q and Q' are small, so that dependence of Cli on ing again the properties of the projection operator from 
either of the bosonic momenta can be neglected. Us- Eq. (|3.56p . we find that 

I 



Tr 



u ur^u U'r^^ ^ S "t 
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16 
1 



rTr 



{(1 + pr) (1 + St) + ip + r){s + t)}I+ {(1 + pr) {s + t) + {1 + st) {p + r)} C 



8 



[{l+pr){l + st) + {p + r) (s + t)] 



(3.62) 
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FIG. 9. The third-order Cooper-channel diagram for the ther- 
modynamic potential. 



This expression vanishes if at least one of the indices from 
the set {p, r, s, t} is different from the others. Therefore, 

only electrons from the same branch contribute to 6^)^ , 
the Zeeman energy can again be absorbed into the chemi- 
cal potential, and 5.=,^ does not depend on the magnetic 
field. 

Finally, the last diagram, 6^e corresponds to the first- 
order self-energy inserted twice into the zeroth order ther- 
modynamic potential. Such an insertion only shifts the 
chemical potential and, for a q dependent U, gives a reg- 
ular correction to the electron effective mass but does not 
produce any nonanalytic behavior. 



IV. COOPER-CHANNEL RENORMALIZATION 

A. General remarks 

The second-order nonanalytic contribution to the spin 
susceptibility comes from "backscattering" processes in 
which two fermions moving in almost opposite directions 
experience almost complete backscattering. Because the 
total momentum of two fermions is small, the backscat- 
tering process is a special case of the interaction in the 
Cooper (particle-particle) channel. Higher order pro- 
cesses in this channel lead to logarithmic renormalization 
of the second-order result i^^— '^^i^^^^^'^^ iSO For a weak in- 
teraction, considered throughout this paper, this is the 
leading higher-order effect. In the absence of the SOI, 
resummation of all orders in the Cooper channel leads to 
the following scaling form of the spin susceptibility 



Sx oc 



E 



In^ (E/A) ' 



(4.1) 



where A is the ultraviolet cutoff and E = maxjT, |A|} 
is small enough so that mU\ ln(_E/A) | ^ 1. As we see, 
both the linear-in-£' term, which occurs already at sec- 
ond order, and its logarithmic renormalization contain 



the same energy scale. The reason for this symmetry is 
very simple: an arbitrary order Cooper diagram for the 
thermodynamic potential contains two bubbles joined by 
a ladder. The spin susceptibility is determined only by 
diagrams with opposite fermion spins in each of the bub- 
bles. Therefore, all Cooper bubbles in such diagrams are 
formed by fermions with opposite spins, so that the log- 
arithmic singularity of the bubble is cut off at the largest 
of the two energy scales, i.e, temperature or Zeeman en- 
ergy. At zero incoming momentum and frequency, the 
Cooper bubble is 



n-^Z^in ^, 

'^ 2n max {T, A} 



rti 



(4.2) 



hence the symmetry of the result with respect to in- 
terchanging T and A follows immediately. It will be 
shown in this section that this symmetry does not hold 
in the presence of the SOI. The reason is that the Rashba 
branches are not the states with definite spins, and di- 
agrams with Cooper bubbles formed by electrons from 
the same branches also contribute to the spin suscepti- 
bility. Although a Cooper bubble formed by electrons 
from branches s and s' 



n^ 



^J2'^ I ^ I d^p9s{^^P)9s'{-^,-p) 



2tt / 2tt 



m , 
— m 



A 



(4.3) 



2tt max{r, |s-s'|Akp} 

looks similar to that in the absence of the SOI [Eq. (|4.2p ] , 
its diagonal element 11^* depends only on T even if T <C 
Akp. Therefore, for A = 0, the Cooper logarithm in Eq. 
(|4.1|) will depend only on T in the limit of T — )• while 
the energy E in the numerator may be given either by T 
or by |a|fci?. 



B. Third-order Cooper channel contribution to Xzz 

Is this section, we obtain the third-order Cooper chan- 
nel contribution to Xzz- This calculation will help 
to understand the general strategy employed later, in 
Sees, live 31 and IIV C 41 in resumming Cooper diagrams 
to all orders. 

The third-order Cooper diagram for the thermody- 
namic potential, depicted in Fig. [HI is given by 



El 

6 



SEil^ = ^ J2 TriGKGpGL) 



K,P,L,Q 



X Tv{G^k+qG-p+qG-l+q)- 

(4.4) 

First, we evaluate the traces Tt::{GkGpGl) — 
J2rst^rstgr{K)gs{P)gt{L) with the coefficients 
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S,,t = Tr[a.(k)n,(p)nt(l)] 



-[AkApAi + irsto^ A(k x p + p x 1 + 1 x k)^ 



4AkApAi 

+ rsia^k. ■ p + A2)Ai + ^(a^k • 1 + A2)Ap + st{aH ■ p + A2)Ak]. (4.5) 



fermionic momenta, 



Since q is small and Brst is an even function of the or the effective Zeeman energy, see Eq. ()4.3p . The third- 
order contribution contains one such logarithmic factor 
which can be extracted from any of the three Cooper 
bubbles; this gives an overall factor of three: 



Tr 



a-'(-k + q)^s' (-P + q)f^f'(-l + q) 



Tr 



Br 



r!w(-k)17,,(-p)a'(-l) 
and the thermodynamic potential becomes 

C/3 



(4.6) 
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EEE 

K,P,Q rst r's'V 

xgr{K)gr'{-K + Q)g,iP)gs>i-P + Q)i4.8) 

V^ 2_. y^ ^rstBr's't'gr{K)gr'{—K + Q) whcre 0ki = ^(k, 1). With this procedure, the third-order 

K,PX,Q rst r's't' diagram reduces effectively to the second-order one, but 

XQ (P)n '(—P + 0)Qt(l)Qt'(—L + O) (4 7) with a new set of coefficients. Since we already know that 

the nonanalytic part of the second-order diagram comes 

Each pair of the Green's functions with opposite mo- from processes with k sa — p, the coefficient Brst (and its 

menta forms a Cooper bubble, which depends logarithmi- primed counterpart) simplify significantly because k x 

cally on the largest of the two energy scales-temperature p = and p x 1 = 1 x k. Integrating over 6'ki, we obtain 
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2tt 
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- 2rstr's't'a*kj,A'^ + -{rtr't' + sts't' - rts't' - r't' st)a'^k%lS.l^ 



A^i^L + (^5 + ^* + ^OA' - rsa^kl][Kl^ + {r' s' + r't' + s't')A'' - r's'a^kl]^ 



(4.9) 



and 



U^ 



5Ei^J = —Y,J2^rstr's't'n*/ J2 9r{K)gr'{-K + Q)g,{P)gA-P + Q)- 



(4.10) 



rst r's't' 



K,P,Q 



The convolution of two Cooper bubbles in Eq. (I4.10p can be re- written via a convolution of two particle-hole bubbles 
by relabeling Q ^ Q + K + P: 

Y, gr{K)gr'{^K + Q)g,{P)gA-P + Q)^Y.Y.9r{K)9s'{K + Q)Y,9s{P)9r'{P + Q)Y.nrs'{Q)n,r'm)\^) 

K,P,Q Q K P Q 

where lists' {Q) is a particle-hole bubble defined in Eq. p.6p . Summing over the Rashba indices, integrating over the 
momentum, and assuming that the Zeeman energy is the smallest energy in the problem, i.e., that A <C max{r, \a\ kp}, 
we find 
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(4.12) 



with J'iy) given by Eq. ([323 • 



The asymptotic behavior of Xzz for T ^ \a\kF is com- puted from Eq. (|4.12p 



^.a>=8x.(^) 



In — 

Ef a 



(4.13) 
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since | ln(T/A)| ^ | ln(|a|fci?/A)|. As it was advertised in 
Sec. IIVAI the TlnT scaling at high temperatures is re- 
placed by the |a| InT scaling at low temperatures which 
implies that the energy scales jaj/ci? and T are not inter- 
changeable. 



FIG. 10. 



Left: the effective scattering amplitude 



riis2;s3S4(k,k'; p, p') in the chiral basis. Right: a 

generic n-th order ladder diagram in the Cooper channel, 



<r,) 



i(k, -k;p, -p) 



C. Re-summation of all diagrams in the Cooper 
channel 



1. Scattering amplitude in the chiral basis 



and, as to be expected, SxzJ scales as TlnT 
In the opposite limit of T ^ |Q;|fci?, 
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(4.14) 



It is more convenient to resum the Cooper ladder dia- 
grams in the chiral basis, in which the Green's functions 
are diagonal. Introducing Rashba spinors |k, s), we re- 
write the number-density operator as 

Pq^Yl H^'^ + "^'*2|k,si)cj^+q^^^Ck,si, (4.15) 

k Si,S2 

so that the Hamiltonian of the four-fermion interaction 
becomes 



^int — 



2 -^ 
q 



Uici)p^p-^ = l J2 J2^i\l;sssAKk';P:P')cl'^ 



p,p',k,k' {si} 



S4'-'P,S3 



Ck,si Ck',S2 : 



(4.16) 



where the effective scattering amplitude is defined by (cf. 
Fig. [TUD 

rW^.,3,^(k,k';p,p') = f/(k-p)(p,S3|k,si)(p',S4|k',S2). 

(4.17) 

In the absence of the magnetic field. 



K 



It.) = ^_ 



where 6*^ = ^ (k,a;) , and Eq. (j4.17p givei 



(4.18) 



.31 



rW^,,3,^(k,k';p,p') = -(7(k-p) 



1 + S2S4eH''p 



1 + 81836 



i(ep-8k) 



(4.19) 



In order to resum the ladder diagrams for the thermo- 
dynamic potential to infinite order, we consider a skele- 
ton diagram depicted in Fig. Ill[ which is obtained from 
the second-order diagram -shown in Fig. [6^~ by replac- 
ing the bare interaction C/(q) with the dressed scatter- 
ing amplitudes: rsjS2;s3S4(k, -k -|- q; p, -p + q) and its 
time reversed counterpart. The dressed amplitudes con- 
tain infinite sums of the Cooper ladder diagrams shown in 
Fig. [ini We will be interested in the limit of vanishingly 
small magnetic fields and temperatures smaller than the 
SOI energy scale: A ^ T <C |a|fci?. In this limit, the 
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FIG. 11. A skeleton diagram for the thermodynamic poten- 
tial H. 




largest contribution to the ladder diagrams comes from 
the internal Cooper bubbles formed by electrons from the 
same Rashba subbands. Each "rung" of this ladder con- 
tributes a large Cooper logarithm L = {m/2TT)\n{A/T), 
which depends only on the temperature, and one has to 
select the diagrams with a maximum number of L factors. 



17 



2. Renormalization group for scattering amplitudes 

Resummation of Cooper diagrams is performed most 
conveniently via the Renormalization Group (RG) 
procedure^. In the Cooper channel, the bare amplitude 
is given by Eq. (|4.19[) with p' = — p and k' = — k or, 
equivalently, ^k' = ^k + "' ^^nd O-pt = 0p + tt: 



riH.;s3s.(k,-k;p,-p) = C/W 



S2;S3S4 



+ V^^,,,,,S0^-<^^ + <^^,3.,e2'(^-«''), (4.20) 

where the three terms correspond to orbital momenta £ = 
0, 1, 2, respectively. The bare values of partial amplitudes 
are given by 



M^.W ;,3,, - (t//4)siS2S3S4. 



(4.21a) 
(4.21b) 
(4.21c) 



Consider now a ladder diagram consisting of n interac- 
tion lines and 2(n — 1) internal fermionic lines, as shown 
in Fig. 1101 As we have already pointed out, in the limit 
T <^ \a\kp, the dominant logarithmic-in-T contribution 
originates from those Cooper bubbles which are formed 
by electrons from the same Rashba branch. Therefore, 
the n-th order Cooper ladder can be written iteratively 
as 

ri;l;.4.3(k,-k;p,-p) = 

= -L f ^ri.'^,-iL(k,-k;i,-i)r«,3,^(i,-i;p,-p) 

(4.22) 

where n > 2. Since only "charge- neutral" terms of the 
type e^^^p^^'^e*'-^'^^'''' survive upon averaging over 9i, dif- 
ferent partial harmonics are renormalized independently 
of each other, i.e., wc have the following group property: 

rll;.3.Jk,-k;p,-p) = (-L)"-Mt^iri;s3s. 
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(4.23) 



Differentiating Eq. (|4.22l) for n = 2 with respect to L we 
obtain three decoupled one-loop RG flow equations 



dL 



UsiS2;S3S4{L) — 7 ^UsiS2-,Ss{L)Uss;S3Si{L), 



(4.24a) 

^-^Vs,s2;S3sAL)^Y.^S4S2-ML)Vss,S3sAL), (4.24b) 



-dZ^' 



SlS2;S3Si 



(L) = yWs,s2;ssiL)W,,,,,sAL) 



(4.24c) 



Solving this system of RG equations and substituting 
the results into the backscattering amplitude. 
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(k, k; k, kj — UsiS2; 
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which is a special case of the Cooper amplitude for p = 
— k, we obtain 



rss;±s±s(k, -k;-k, k) 
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r,(±^.^^)(k,-k;-k,k)=0. 



(4.26) 



(4.27) 



(4.28) 



where cr(siS2; S3S4) = {(S1S2; S3S4), (s2, S3; S4S1), 
(s3, S4; siS2)j (s4, Si; S2, S3)} stauds for all cyclic permu- 
tations of indices. 

We see that the RG flow described by Eqs. (|4.24at - 
I4.24cp has a non-trivial solution: whereas the amplitudes 
in Eqs. (|4.26p and (|4.28p flow to zero in the limit of 
L — > cx), the amplitudes in Eq. (|4.27p approach RG- 
invariant values of ±C//2. This behavior is in a striking 
contrast to what one finds in the absence of the SOI, 
when the repulsive interaction is renormalized to zero in 
the Cooper channel. Notice that we consider only the en- 
ergy scales below the SOI energy, while the conventional 
behavior is recovered at energies above the SOI scale. 

The scattering amplitudes can be also derived by iter- 
ating Eq. (|4.22p directly. Examining a few first orders, 
one recognizes the pattern for the n-th order partial am- 
plitudes to be 
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(4.29) 
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rl(UTT)(k.-k;-k,k).= o. 



(4.30) 
(4.31) 



Summing these amplitudes over n, one reproduces the 
RG result. 



3. Renormalization of Xzz 

The infinite-order result for the thermodynamic po- 
tential is obtained by replacing the bare contact interac- 
tion U by its "dressed" counterpart F in the second-order 
skeleton diagram Fig. [11] 



with initial conditions specified by Us-^s2:s3S4,iO) = '^"^x — ^t 2^ 2^rsjs^;s3s.j(k, — k; k, k) 



f^SlS2;S3S4J '''^SiS2;S3S4 (L") — ''S1S2;S3S4 7 and \/sj^s2;S3S4 (Uj — 

yil) 

l'siS2;S3S4 ■ 



^ f s3S2;siS4(^k, k; k, — k)IIsjS2f-ls3S4 (4-32) 
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with r 



S3Si:SiS2 



(-k,k;k,-k) 



SlS2;S3S4, 



(k,-k;-k,k). 



Now, we derive the asymptotic form of 6xzz vahd in the 
hmit of strong Cooper renormalizaton, i.e., for UL 3> 1. 
In this hmit, the only non- vanishing scattering amplitude 
is given by Eq. (g^. Replacing the full T by its RG- 
invariant asymptotic limit rs_s;ips±s(k, — k; — k, k) = 
±C//2, we obtain 



^-.^-^ri:/$[,n5 
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(4.33) 
In contrast to the perturbation theory, where the 
magnetic-field dependence of the thermodynamic poten- 
tial was provided by the vertices while the polarization 
bubbles supplied the dependence on the temperature and 
on the SOI, the vertices in the non-perturbative result 
(|4.33p depend neither on the field nor on the SOI. There- 
fore, the dependences of 2 on all three parameters {B, 
T, and a) must come from the polarization bubbles. The 
integral over q along with the sum over the Matsubara 
frequency ft have already been performed in Sec. IIIII 
Note that the last term in square brackets (proportional 
to Ilg) does not depend on the magnetic field and thus 

I 



can be dropped. The final result reads 
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where use was made of the expansion 
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(4.34) 



(4.35) 



(4.36) 



and made use of the asymptotic form (I3.27P of the func- 
tion F to find that T'{x) ~ x^ for a; ^ 1. Comparing the 
non-perturbative and second-order results for Xzz [given 
by Eqs. (|4.35p and p.30p . respectively], we see that the 
only effect of Cooper renormalization is a change in the 
numerical coefficient of the nonanalytic part of Xzz ■ This 
is a consequence of a non-trivial fixed point in the Cooper 
channel which corresponds to finite rather than vanishing 
Coulomb repulsion. 

The temperature dependence of Xzz can be also found 
for an arbitrary value of the Cooper renormalization pa- 
rameter UL. Deferring the details to Appendix IC 1[ we 
present here only the final result 
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(4.37) 



In the limit of strong renormalization in the Cooper chan- 
nel, i.e., for UL ^ 1, only the first term survives the loga- 
rithmic suppression, and Eq. (|4.37p reduces to Eq. (|4.35p . 
In the absence of Cooper renormalization, i.e., for L = 0, 
Eq. (|4.37p reduces to the second-order result p.30p . In 
between these two limits, 5xzz is a non-nionotonic func- 
tion of UL: as shown in Fig. [121 ^Xzz exhibits a (shal- 
low) minimum at UL w 2.1. In a wide interval of UL 
(0.9 < UL < 5.6), the sign of Sxzz is opposite (negative) 
to that in either of the high- and low-temperature limits. 
It is also seen from this plot that the low-T asymptotic 
value (marked by a straight line) is reached only at very 
large (> 100) values of UL. 



4-. Renormalization of Xxx 

The second-order result for the in-plane magnetic 
field is renormalized in a similar way with two excep- 
tions. First, because the Zeeman energy is anisotropic 



I 

in this case - the effective magnetic field A±]f^p{d\) = 
\J a^k\ ± 2afcFAsin6'i -I- A^ depends on the direction 
of the electron momentum 1 with respect to the field 
Q\ = Z(l, B)- integration in the "rungs" of the Cooper 
ladder can be performed only over the fermionic fre- 
quency and the magnitude of the electron momenta (or, 
equivalently, the variable ^i). Consequently, the elemen- 
tary building block of the ladder 



o-^E 



ra 

2^ 



m 

— -—In 

27r max 



c?65t(w,l)gt'(-t^,-l) 



A 



{T,|tAk^(0i)-t'A_k^(0i)|} 

(4.38) 

depends on 6i. 

In principle, the dependence of L on the angle 9\ 
should be taken into account when averaging over 0i. 
However, in the limit of A <C T <C |q:|A:f, the angle- 
dependent term under the logarithm can be approxi- 
mated as \tAkp{6i) — t'A_fc^(0i)| w |t — t'l \a\kF and, as 



19 




FIG. 12. Nonanalytic part of x^^^^-: normalized by the second- 
order result (|3.30|) . as a function of the Cooper-channel renor- 
malization parameter UL = {Um/2TT)\n{A/T). The horizon- 
tal line marks the low-temperature limit. 



it was also the case for Xzz, L — (m/27r) ln(A/T) pro- 
vided that t = t' . 

Second, the particle-hole bubbles also depend on the 
direction of the electron momentum, hence, the infinite- 
order result for the thermodynamic potential is found 
again by replacing the bare interaction U in the second- 
order diagram by rsjS2;s3s^ {9^} and retaining the angular 
dependence of the bubbles. 

The RG equations for the in-plane magnetic field are 
considerably more complicated. The main difference is 
that even the RG-invariant terms depend on the mag- 
netic field. A detailed discussion of Cooper renormaliza- 
tion of the scattering amplitudes and spin susceptibility 
for this case is given in Appendices IC 2al and IC2b[ re- 
spectively. Below, we only show the final result for the 
renormalized spin susceptibility 



Xo fmll 



Ifcf 



Ep 



T 



O TZT^ ■ (4-39) 



The T-independent term is the same as without the 
Cooper renormalization [cf. p.5ip ]. The linear- in-T term 
however is suppressed by at least a factor of 1/ In T, sim- 
ilar to the case of no SOI, where it is suppressed by a 
In^ T factor. 



V. SUMMARY AND DISCUSSION 

We have considered a two-dimensional electron liq- 
uid in the presence of the Rashba spin-orbit interaction 
(SOI). The main result of this paper is that the combined 
effect of the electro-electron and spin-orbit interactions 
breaks isotropy of the spin response, whereas either of 
these two mechanisms does not. Namely, nonanalytic 
behavior of the spin susceptibility, as manifested by its 
temperature- and magnetic-field dependences, studied in 
this paper, is different for different components of the sus- 
ceptibility tensor: whereas the nonanalytic behavior of 



Xzz is cut off at the energy scale associated with the SOI 
(given by \a\kF for the Rashba SOI), that of Xxx (and 
Xvv = Xxx) continues through the SOI energy scale. The 
reason for this difference is the dependence of the SOT 
induced magnetic field on the electron momentum. If 
the external magnetic field is perpendicular to the plane 
of motion, its effect is simply dual to that of the SOT 
field: the T dependence of Xzz is cut by whichever of 
the two fields is larger. If the external field is in the 
plane of motion, it is always possible to form a virtual 
particle-hole pair, which mediates the long-range inter- 
action between quasiparticles, from the states belonging 
to the same Rashba branch. The energy of such a pair 
depends on the external but not on the effective field, so 
that the SOI effectively drops out of the result. We have 
also studied a non-perturbative renormalization of the 
spin susceptibility in the Cooper channel of the electron- 
electron interaction. It turns out the RG flow of scat- 
tering amplitudes is highly non-trivial. As a result, the 
spin susceptibility exhibits a non-monotonic dependence 
on the Cooper-channel renormalization parameter (In T) 
and eventually saturates as a temperature-independent 
value, proportional to the SOI coupling \a\. 

Notably, all the results of the paper are readily ap- 
plicable to the systems with large Dresselhaus SOI and 
negligible Rashba SOI. In this case the Rashba spin-orbit 
coupling should be simply replaced by the Dresselhaus 
spin-orbit coupling. 

Now we would like to discuss possible implications of 
these results for (in)stability of a second-order ferromag- 
netic quantum critical point (QCP). This phenomenon 
depends crucially on the sign of the nonanalytic correc- 
tion. In this regard, we should point out that we limited 
our analysis to the simplest possible model, which does 
not involve the Kohn-Luttinger superconducting instabil- 
ity and higher-order processes in the particle-hole chan- 
nel. Therefore, the sign of our nonanalytic correction is 
"anomalous", i.e., the spin susceptibility increases with 
the corresponding energy scale. As in the absence of the 
SOI, however, either of these two effects (Kohn-Luttinger 
and particle-hole) can reverse the sign of nonanalyticity. 
Therefore, it is instructive to consider consequences of 
both signs. 

In the absence of the SOI, a nonanalyticty of the 
anomalous sign renders a second-order ferromagnetic 
QCP unstable with respect to either a first-order phase 
transition or a transition into a spiral state 4i^ This result 
was previously believed to be relevant only to systems 
with a SU{2) symmetry of electron spins; in particular it 
was shown in Ref. i5ib that there is no nonanalyticity in x 
for a model case of the Ising-like exchange interaction be- 
tween electrons. We have shown here that broken (by the 
SOI) SU{2) symmetry is not sufficient for eliminating a 
nonanalyticity in the in-plane component of the spin sus- 
ceptibility (xxx)- Based on our results for the magnetic- 
field dependence of Xzz and Xxx [cf. Eqs. p. 371) and 
(I3.53|) ] , we can construct a model form for the free energy 
as a function of the magnetization M . The most inter- 
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esting case for us is the one in which the Zeeman energy 
due to spontaneous magnetization, ~ M/m^B is larger 
than the SOI energy scale |a|/c_F, so that Xxx 7^ Xzz- Ig- 
noring Cooper-channel renormalization, we can write the 
free energy as 



aM^ -b{\Ma, 



\M„ 



M^ 



(5.1) 



.1/2 



where M = (M^ + M^ + M^) "' " and the coefficient of 
the quartic term was absorbed into the overall scale of 
F^ which is irrelevant for our discussion. An important 
difference of this free energy, compared to the case of 
no SOI, is easy-plane anisotropy of the nonanalytic, cu- 
bic term. In the absence of the cubic term (6 = 0), 
a second-order quantum phase transitions occurs when 
a = 0; in the paramagnetic phase, a is positive but small 
near the QCP. Since x is isotropic at the mean-field level 
[cf. Eq. ([2^ ]. the regular, M^ and Af* terms in Eq. ([5J|) 
are isotropic as well. If 6 > (which corresponds to the 
anomalous sign of nonanalyticity) , the cubic term leads 
to a minimum of F at finite M; when the minimum value 
of F reaches zero, the states with zero and finite magne- 
tization become degenerate, and a first-order phase tran- 
sition occurs. The first-order critical point is specified by 
the following equations 



dF 
dlVL 



0, 



dF 



0, 



dF 
dMy 



0, and F = 0. (5.2) 



For a > 0, the only root of the first equation is Mz = 0, 
i.e., there is no net magnetization in the z direction. 
Substituting Mz = into the remaining equations and 

I 



employing in-plane symmetry (Mx = My), we find that 
the first-order phase transition occurs at a = 6^/8. The 
broken-symmetry state is an XY ferromagnet with spon- 
taneous in-plane magnetization Mx — My' — 6/4. 
The first-order transition to an XY ferromagnet occurs 
if the Zeeman energy, corresponding to a jump of the 
magnetization at the critical point, is larger than the SOI 

energy, i.e., Mx /mjiB ^ |a|fcF- In the opposite case, 
the SOI is irrelevant, and the first-order transition is to 
a Heisenberg ferromagnet. 

If the nonanalyticity is of the "normal" sign {h < 0), 
the transition remains second order and occurs at a = 0. 
However, the critical indices are different for the in-plane 
and transverse magnetization: in the broken-symmetry 
phase (a < 0), Mx = My ex (—a) while Mz ex (— a)^/^. 
Since \a\ ^ 1, the resulting state is an Ising-like ferro- 
magnet with Mz^ Mx = My. 

A detailed study of the |q| dependence of x in the 
presence of the SOI will be presented elsewhere.— 
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Appendix A: Temperature dependence of the spin susceptibility of free Rashba fermions 

In this Appendix, we consider the temperature dependence of the spin susceptibility of free 2D electrons in the 
presence of the Rashba SOI. The transverse and parallel spin susceptibilities, x^z a-nd x^x , are given by 



X 



{0} _ 



-Y.Ty[g{K)&zG{K')&z 

K 
K 



(Al) 
(A2) 



where K = (w, k) and K' = (w, k -f q) with g — >► 0. Evaluating the traces, we obtain for x 
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.{0} 



--^E/7S2E^(l-^^)5«(-'k)5*(^,k)^-2rE|^.9+(-,k)g_(c.,k), 



(A3) 



where we took advantage of the isotropy of g±(w, k) and put q = Q, because the poles in the Green's functions of 
different branches reside on opposite sides of the real axis. We see that Xzz is determined only by inter-subband 



transitions. Similarly, we obtain for Xxx 
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^J2 [ 7^ E ^ [1 - sicos(20k)]<?.(w, k)gtiuj, k + q)|,^o. 
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Since cos6'k averages to zero, Xxx can be written as 



.{0} ^ 1(0} , ij {0} 



(A5) 



where 



/72 7 
— -2[5+(w,k).9+(c^,k + q)+.g_(a.,k).9_(a.,k + q)]|g^o 
(ztt) 



{0} 



is the contribution from intra-subband transitions, absent in x 

Let us evaluate Xzz first. Performing the fermionic Matsubara sum, we obtain 

1 1 



^{0} 

A,zz 



f°° kdk v-^ 1 

Jo 2n ^^ 2ak 



i^ ^ £,k ~ otk icj — (^k + Oik 



(A6) 



1 f°° 

/ dk[nF{^k~ak)~nF{^k + ak)] (A7) 

27ra Jo 



with a Fermi function npie) — [e*-'^^'/^ + l] and ^k — k^/2m — fi. Changing variables in the first integral to 
^k ^ ckfc = e_ — II, we find two roots: k_ = ma — ^J{raa)'^ + 2me^, valid for — eq < e_ < with dk_ /de- < 0, 
and fc_ — ma + y^map + 2meZ, valid for e_ > with dk_ /de- > 0, where eq = m.a^/2. Similarly, we change 
variables in the second integral to ^k + ctk = e+ — /i and obtain only one positive root fc+ = —ma + y/{ma)'^ + 2me^, 
valid for e+ > with dk+/de+ > 0. Notice that the absolute values of the (inverse) group velocities are the same for 
both branches: |dfcl' /de_| = |dfc+/de+| = m[{ma)'^ + 2me±]^^/'^. Therefore, 
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(A8) 



where we dropped the index on the integration variable e. Notably, the high energy contributions from the two Rashba 
branches cancel each other for any value of the chemical potential and the spin susceptibility is determined exclusively 
by the bottom part of the lower Rashba branch. Integration by parts yields 



d"J^Xo[nF{0)- 1^ de^ 



e dnp[e) 
eo de 



(A9) 



where 71^(0) = [e~'^/'^ + 1]~^. In order to evaluate this integral, it is convenient to consider three limiting cases. 

For T <C eo ^ /i, i.e., when both Rashba subbands are occupied and the temperature is lower than the minimum 
of the lower subband, we approximate rii?(0) w 1 — e^^^^, —TdnF{e)/de « e^'^^^/^, and ^JT+TJeQ sa 1 + e/2eo, so 
that 
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At T = 0, xl°^ = Xo. ^ 

For eo <^ T <^ fi, i.e., when again both Rashba subbands are occupied but the temperature is higher than the 
minimum of the lower subband, we keep the same approximations for np(0) and —Tdnp{e)/ de but neglect the e 
dependence of the Fermi function in the integrand: 
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(A12) 



Finally, for /i < 0, i.e., when only the lower Rashba subband is occupied, the first term in Eq. (|A9[) gives only 
an exponentially weak temperature dependence, while the Sommerfeld expansion of the second term generates a T^ 
contribution because the density of states depends on e: 
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(A13) 
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At zero temperature, Xzz = Xo\/^ ~ ImI/^o vanishes at the bottom of the lower subband. 
We now calculate Sx%J given by Eq. (jA6p . which can be written as 
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The change of variables is straightforward in the second integral, since the equation e+ — /i = ^k + oik has only one 
positive root fc+ and the density of states is given by 
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with eo = vno? j2. In the first integral more care must be taken: for e > 0, we have k_ = men + \J {ma)'^ + 2me and 
the density of states is given by 
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however, for — eo < e < both roots enter the density of states 
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Summing up the contributions from all energies, we find 
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where the second and third integrals are easily evaluated because e drops out of the sum v^-\- v^ ~ tr/tt. Combining 
the above result with Eqs. (IA2p and (|A8p . we get 
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For r < eo < /i, we approximate 77f(0) w 1 - e"''/^ and -~TdnF{e)/de = e(e-A')/T(g(£-M)/T _^ ^)-2 _ g(e-M)/T ^ 
before, and expand (1 + e/2eo)/^l + e/eg ~ 1 + e^/Seg, so that 
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For eo <C T ^ /i, we keep the same approximations for the Fermi function but, as it was also the case for Xzz , neglect 
the e dependence of the Fermi function in the integrand 

„ /" ,,l±i^^Mf) «le-./- rde^^^BLe-^'-^le-^l- H de^^^^L = '-^e-^'- . (A21) 
i-.n ^^+7l^. de T Jo ^/T^TJTo T J^ ^/^^77^ 3T 
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Finally, for // < 0, the Sommerfeld expansion of the second term in Eq. (IA19P yields 
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for T <ti min{|^|, Co — |^|}. At zero temperature, Xxx = Xo (1 — lMl/2eo) / (1 — ImI/^o) diverges at the bottom of 
the lower subband. 
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Appendix B: Absence of a small-g singularity in the static polarization bubble with spin-orbit interaction 

In Sees, mil and fill D[ we argued that there is no contribution to the nonanalytic behavior of the spin susceptibihty 
from the region of small bosonic momenta, q <^ kp- This statement contradicts Ref. [32 where it was argued that, in 
the presence of the SOI, a static particle-hole bubble has a square-root singularity a.t q = qq = 2ma (in addition to the 
Kohn anomaly which is also modified by the SOI). For a weak SOI, go is much smaller than kp and thus the region 
of small q may also contribute to the nonanalytic behavior. Later on, however, Refs. '28! andlSJ showed that there is 
no singularity a,t q = qo- According to Ref. [2^, the reason is related to a subtlety in approaching the static limit of a 
dynamic bubble. While we agree with the authors of Refs. [2^ in that there are no small- g singularities in the bubble, 
we find that a cancellation of singular terms occurs in the calculation of a purely static bubble. The same result was 
obtained in an unpublished work.'^* For the sake of completeness, we present our derivation in this Appendix. 

Evaluating the spin trace, we obtain for the static polarization bubble 

n(g) = Yl TriG^^.k+qGc^.k] = 2 X!X! [-^ + st cos((/3k+q - -^k)] 9s{i^, k + q)gt{uj, k), (Bl) 

K K s,t 

where, as before, K — (a;,k), v'k-i-q = ^(k + q, e^,) and (/Sk = ^(k, e^,). We divide n(g) into intra- and intersubband 
contributions as 



where 



n(<z) = n++(g) + n__(g) + n±(g), (B2) 

n±±(g) = :t X^t^ + cos((^k+q - tpk)]g±iuj, k + q)g±(w, k), (B3a) 

^ K 

n± (g) = - ^ [1 - cos(v3k+q - fk)] [g+ (w, k + q)5_ (w, k) + g_ (w, k + q)g+ (w, k)] 
^ K 

= X][l ~ cos((^k+q - (pk)].g+(w, k -I- q)g-(a;, k). (B3b) 

K 

In the last line, we employed obvious symmetries of the Green's function. 

First, we focus on the intersubband part, Ti±{q). Summation over the Matsubara frequency yields 

T2_^g+{p,uj)g^{k,uj) = ^ — — , (B4) 

^P ^k 

where p = k + q, ^^ = ^i^ ± ak and, as before, ^k = k^/2m — Ep. Introducing additional integration over the 
momentum p, as it was done in Ref. 32, Eq. (jB3bp can be re-written as 
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where 6* = Z(k, q). Integration over 9 yields 



/ d0S{p^ -k^ -q^ - 2kq cos e)(l- 



kq cos 9\ 1 q'^ — {k — pY 



kp J kp ^(k + g)2 - p2^p2 - (fc - g)2 ' 



(B6) 



which imposes a constraint on the range of integration over p, i.e., |/c — g| < p < k + q. Since we assume that 
q <^ k ^ p~ kp, Eq. (|B6p can be simplified to 
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For a weak SOI (m|a| ^ kp) , £^^ sa ^k i akp. Switching from integration over k and p to integration over ^k and ^p, 
we find 

n± g) ^ ,23, / rfCk / c?Cp J(wFgr- (Ck-Cpr ? — . , r, i. • (^9) 

An-^v'pkF J .^ J^^-v^q ^ 4p-4k + 2afcj. 

Shifting the integration variables as Cp ~^ ?p ~ afc^, Ck ^^^ Ck + otkp, we ehminate the dependence of the Fermi 
functions on akp. Assuming also that T = 0, we obtain 

n±(9) = , 2 3. / '^^k / dCp v/(«f 9)' - (^P - ek - 2akFY ""^ ^P^ ^^ ^'^^ (BIO) 

'ill VpKp J ^oo J(,i^+2akF-vpq ^ ?p — ?k 

where Q (x) is the step function. Notice that the integrand is finite only if ^kCp < Oi which imposes further constraints 
on the integration range. 

We will now prove that ^±((l) given by Eq. (|B10[) is continuous at g = go- To this end, it is convenient to consider 
the cases of g < go and q > qq. Combining all the constraints together, we find that Il± for q < qq can be written as 
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VFq—2akF J ^k.-\-2akF—VFq J >P "^k 

Reversing the sign of ^k, absorbing ^k into ^p, and defining the dimensionless variables x — ^-^jvpq and y — fp/wF?, 
we obtain 

where /3 = qo/q > 1. Next, we switch the order of integration in the first term, so that the integrals over x can be 
readily evaluated 
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For q > qo, we have 
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Manipulations similar to those for the previous case yield 

47r^i;FfcF ^0 Jx y Jo Jx y J 

with /? = qo/q < 1. Interchanging the order of integrations over x and y, we find 

n,(g>2^«) = -^^"^,yr3(^T^ + £"^,yr3(^3^^ _ (bi6) 

Since n^(g = go — 0+) = ^±{q — qo + 0^), the function n±(g) = —{m/2Tr){q/2kp)'^ is continuous at g = go and, 
thus, there is no singularity in the static particle-hole response function. In addition, n-|-(g) does not depend on the 
SOI. However, since there is no g^ term in the 2D bubble for g < 2kp in the absence of the SOI, the g^ term must be 
canceled out by similar terms in 11++ (g) and 11 (g), which is what we will show below. 
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Having proven that n±(g) is an analytic function of q, we can re-derive its q dependence simply by expanding 
the combination (pk+q — ^'ia in Eq. (IB3b[) for q <^ kp &s. ^p\^+c^ — Vk ~ (^z/^f) siniyj^qj where (/Jkq = ^(k, q); then 
1 — cos(i^k+q — <(2k) ~ {q/kpY sin ((ySkq) /2. Sincc we already have a factor of q^ up front, the Green's functions in 
Eq. (jB3bp can be evaluated at <? = 0. Accordingly, Eq. (jB3bp becomes 



n±(g) 



2 Vfc 



^ sin^ <y5kq.9+ (k, a;)g_ (k,c 



m ( q \ f npj^j,) - UFJ^^,) 

^ l^J J '^"^ C-Ck 






m I q 
2tt \2kF 



£ u^ ■ (B17) 



Expanding Eq. (jB3a|) for the intraband contribution to the bubble also to second order in q, we obtain 
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27r 



■Sk+q "Sk 



q^O 






Thereby the total bubble ([62 



is independent of q for q <C 2/0^-. 



n(g) 



m 

IT 



(B18) 
(B19) 

(B20) 



Appendix C: Renormalizaton of the scattering amplitudes in the Cooper channel in the presence of the 

magnetic field 

In this Appendix, we present the derivation of the RG flow equations for the scattering amplitudes in the Cooper 
channel in the presence of the magnetic field. These amplitudes are then used to find the non-perturbative results for 
the spin susceptibility. Since the amplitudes are renormalized quite differently for the field applied perpendicularly 
and parallel to the 2DEG plane, we will treat these two cases separately. 



1. Transverse magnetic field 

a. RG flow of the scattering amplitudes 
If the magnetic field is transverse to the 2DEG plane, B ~ Be^, the eigenvectors of the Hamiltonian (|l.ll) read 

1 / (A - sAk)ie-*''Vafc 



|k,s) 



V^Uk) 



1 



(CI) 



where Ns{k) = 2 + 2A(A — sAk)/(afc)^ is the normalization factor and, as before, Ak = (A^ + a^k^)^^/^ is the 
effective Zeeman energy. Since, by assumption, jaj/ci? <^ Ep, we approximate Ak by Akp. Substituting the above 
eigenvectors into Eq. (|4.17p . we find the scattering amplitude 



/ 



r«^^,^,3(k,k';p,p')-C/ 



\ 



n 



'■=1.2 



2A(A-...iAkj,: 



a2fc| 



-(A-siAkJ(A-S3AkJe^(^-- 



(A-s2Ak,)(A-s4A,Je*^V-''.' 



i(V-V) 



(C2) 
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To find the spin susceptibility, we need to know the scattering amplitude to second order in the magnetic field. 
Expanding Eq. (|C2p to second order in S = A/|a|A:F (note that s^ = 1 and s^ = Si) and projecting the amplitude onto 
the Cooper channel, where the momenta are correlated in such a way that k' — — k and p' = — p or, equivalently, 
^k' = ^k + "■ and 6p' = 6p + tt, we obtain 

f(i) ('k -k- n -n^l — r/ A- V p*(^p-^k) , w 2i(ep-ek) (f<n\ 

where we introduced partial amplitudes 

UsiS2;s3Si = y + -^-(si + S2 + S3 + S4)S + T^(siS2 + S1S3 + S1S4 + S2S3 + S2S4 + S3S4 - 2)(5^ + 0{6^), (C4) 

l/siS2;S3S4 =-r(SlS3 + S2S4J + — (S1S2S3 + S1S2S4 + S1S3S4 + S2S3S4 - Si - S2 - S3 - S4)0 
4 o 

+ — (1 - S1S2 - S2S3 - S3S4 - S1S4 - S1S3 - S2S4 + siS2S3S4)^^ + 0{S^), (C5) 

.T. U U, 

WsiS2;s3Si =-rSlS2S3S4 - — (S1S2S3 + S1S2S4 + S1S3S4 + 328354)0 
4 o 

+ T^('Sl'*2 + S1S3 + S1S4 + S2S3 + S2S4 + S3S4 - 2siS2S3S4)5^ + 0{S^). (C6) 

For (5 = 0, the partial amplitudes reduce back to Eqs. (|4.21all4.21cl) . The RG flow equations for the partial amplitudes 
are the same as in the absence of the magnetic field and are given by Eqs. (|4.24al) ~ (|4.24cp with initial conditions 
(|C4p - (|C6[) . Since the differential equations for U, V, and W are identical, for the sake of argument we select the first 
one, copied below for the reader's convenience, 

~ "7T^S1S2;S3S4(-^) = /_^UsiS2;Ss{L)Uss;.S3S4{L) (C7) 

and introduce the following ansatz 

'JsiS2;s3S4[-t^) - l + bL ' 

which satisfies the initial condition. Substituting this formula into the differential equation for U and multiplying the 
result by (1 + bL)^, we obtain an algebraic equation for as-i^s2\s3S4 and b 

^f^SiS2;S3S4(0) - as,S2;S3S4 + ^{UsiS2:Ss{0) + as,S2:SsL){Uss;S3S4{0) + ass-S3S4L) = 0. (C9) 

S 

Grouping coefficients of a polynomial in L, we obtain the following set of equations 

6f/,,,,;,3,,(0) - a.,.2;.3.4 + Y. Us,S2;SsmUs.-S3S4{0) = 0, (CIO) 

s 

Y[Us^S2;ss{0)ass-s3S4 + asiS2;ssUss:S3S4{0)] = 0, (Cll) 

s 

/ ^ ClsiS2.SsO'SS-,S3S4 =0, (C12) 

s 

which are not independent. Thereby, we choose two out of three equations, namely, Eq. (JCIOP and Eq. (IC12[) with the 
si = S2 = S3 = S4 = 1 combination of the Rashba indices, so that there are 17 equations for 17 unknown variables: 
as-^s2;s3S4 and b. The final solutions are listed below 

c/±±;±±(i) = 7i + (1 + 52)^^^/2' (C13) 

U.s;-s-s{L) = Us-s;-ss{L) = Us-s;s-s{L) = ^^-^1-^^^^^, (C14) 

t/.(±T;TT)W - Ti + i^i + 5^)UL/2- ^^^^^ 
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The same procedure is repeated to derive the V- and VF-amphtudes Hsted below 






and 






Summing up aU the contributions to the backscattering amphtude, obtained from the Cooper ampUtude for a special 
choice of the momenta p — — k, i.e., for 0p — 0k = tt, we obtain 



Ts,s2:s,sAK-'i^-,-K'ii) ^ U,,s,;sssAL) -Vs,s2:s,sAL) +Ws,s,-s,sAL). (C24) 

To second order in the field, the backscattering amplitudes read 



r.,.(k, -k; -k, k) ^ (^ - ^^^) + (^(^^ + (^^) ^^ (C25a) 

r.,....(k, -k; -k, k) = (^ + ^^) - (^^^^ - ^^^ + ^) .^ (C25b) 

r__.(k,-k;-k,k) ^ (^ + ^) - (^^ - ^^^ + ^) .^ (C25C) 

r.„_(k,-k;-k,k) ^ (^ - ^) - (^^ - ^^ + ^) .^ (C25d) 

r.(,,^,,, (k, -k; -k, k) ^ T (^^^ + ^) ^. (C25C) 



For 5 = 0, we reproduce Eqs. (J4.26II4.25)) . In the limit of strong renormalization in the Cooper channel, i.e., for 
UL — > oo, the field-dependent terms in Eqs. (jC25allC25eP vanish. A posteriori, this explains why we could obtain the 
renormalized result (j4.35p for Xzz in the main text using only the zero-field amplitudes. 
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b. Renormalization of Xzz in the Cooper channel 



As in Sec, live 3[ the thermodynaniic potential in the presence of Cooper renormahzation is obtained by substituting 
the renormahzed scattering amphtudes (jC25allC25ep into Eq. (|4.32l) 



^S, 



¥^ 



A2 



a^kl ^ 



qdq 

1-K 



1 



u 



u 



2 \2 + UL 



(n^ 



n'-+ - 2n2) 



u 



u 



2 + UL 2(1 + UL) 



(n+_n_+-n2) 



(C26a) 



C/2(16 + 32UL + 22C/2L2 + 6U^L^ + W^L 



4t4:\ 



^^ r-B 



u 



4(l + t/L)2(2 + C/L)2 

2 

^o(^^ 



-n 



u 



2{1 + UL) 2 + UL 
I f U 2U 



H- 



2U 



2no) 
u u 



2 \2{1 + UL) (2 + C/L)2 2 + ULj \2 + UL 2 

U 2U 2U \ f U U 



(n 



tf + - 2n2) 



(C26b) 



2(l + t/L)2 (2 + C/L)2 2 + ULJ \2 + UL 2{l + UL) 
U'^L^{12 + 18UL + 7C/2L2) ^^ 
2(H-C/L)2(2 + t/i)= 



(n+_n_+-n2) 



-nt 



(C26c) 



The first part of S'E.^z (|C26aP came from the field-independent terms in the scattering amplitudes. This part depends 
on the magnetic field through the combinations of the polarization bubbles. The second part (|C26c[) already contains 
a field-dependent prefactor (A2) resulting from the field-dependent terms in the scattering amplitudes. Therefore, 
the polarization bubbles in this part can be evaluated in zero field. The integrals over q along with the summation 
over the Matsubara frequency H. have already been performed in Sec. IIIII Note that in each square bracket the last 
term (proportional to IIq) depends neither on the field nor on the SOI. In fact, it can be shown^i^^ that this formally 
divergent contribution has a cubic dependence on temperature, TJ2nI l^l^o ^ ^^i thus it adds a higher order 
correction and can be dropped. The final result reads 



fczz = 




(C27) 



The spin susceptibility is obtained by expanding Eq. (|C27|) further to order A2. We also need to recall that our 
treatment of Cooper renormalization is only valid for T <C lajfci?, because we kept only the T- but not a-dependent 
Cooper logarithms (see the discussion at the end of Sec lIVCTj) . Therefore, the function T and its derivative should 
be replaced by their large-argument forms, Eq. p.27p . Doing so, we obtain the final result for the nonanalytic part 
of the Xzz, presented in Eq. ()4.37p in the main text. 



2. In-plane magnetic field 

a. RG flow of the scattering amplitudes 



For the in-plane magnetic field, B = Be^, the RG fiow of the scattering amplitudes is more cumbersome because 
the eigenvectors of Hamiltonian (II. ip depend in a complicated way on the angle between the magnetic field and the 
electron momentum 



|k,s) = J-f^^kAA-ze'^'-afc) 



(C28) 
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where Ak = (A^ + 2afcAsin6'k + a^fc^)"^/^ is the effective Zeeman energy. For that reason, the (double) Fourier 
series of the scattering amplitude (|4.17p in the angles ^k and 9p contains infinitely many harmonics: 

oo 

rW.;.3s.(k,-k;p,-p)= Y. ^«4"3"le™'''e"'^ (C29) 



m,n— — oo 



To second order in the field, however, the number of nonvanishing harmonics is limited to 15. Indeed, expanding the 
eigenvector (|C28I) to second order in <5 = A/\a\kp as 

\]^ g\ J_ ( sie-^'^'- [1 - i cos O^S - (cos^ 6'k - i sin 26'k)(5^/2] A _j_ Q/gS-s rQ^Q\ 



V2\ 
and substituting Eq. (jC30[) into the scattering amplitude (|4.17p with k' = — k and p' = — p, we obtain 



r«.;.3.4(k,-k;p,-p) = Us,sr..s,s,+Vs,sr,sss,e''^e-''- + W, 



g2i6lpg-2iek 



Here, 



SiS2;S3S4 

A piO-PpiS-k. I f> p-iSpp-iBk i p 2i9p -iO^^ , ^ iflp -2jek 

SiS2;S3S4^ ^ T^ ^-*SiS2:SsS4^ ^ ^ -'- SiS2;S^S4^ ^ T^ ^^Si S2 ;S3 S4 '^ ^ 

7 p-«^k I T p2iep , 71^ p-2jek 

p p3i9p -iSk I /Q piSpp-Si^k < f? p4iep -2i0k I o p2i0p -4i0k i /-Of/rSx (-rjQ-i ^ 

SiS2;S3S4C '^ T^ ^:siS2 ;S3S4'^ ^ T^ -'"'SiS2;S3S4'^ C T^ ^SiS2;S3S4'^ c n^ \^\^U y. ^^VyO-Ly 



f^s4S2;.3.4 = ^i'M%% = iU/16)[A + (siS3 + S2S,)S% (C32) 

V;4.2;a3a4 = ^i'M'^s-^s} = (f//8)(siS3 + S2S4)(2 - ,52), (C33) 

T^.i.2;.3.4 = rW^^il^ = (C//16)[4siS2S3S4 + (S1S3 + S2S4)^2], (C34) 

^.i.2;s3.4 = r(iHi^i}^ = (C//32)(siS3 + S2S4)S^ (C35) 

R — p(l){-l.-l} — A ('CSfi'l 

^SiS2;S3S4 ~ -*- SiS2;S3S4 — ^SiS2;S3S4, \VjO\JJ 

Fs,S2;s,S4 = rW^^.-^i; - (t//8)i(siS3 - S254)5, (C37) 

(7 _ p(i){i,-2} _ _p CCSS'l 

^SlS2;S3S4 ~ -■- SiS2;S3S4 ~ -f^ Si S2 ;S3S4 : ^_V.yOO; 

ff — p(i){i,o} — _p ('csq'i 

^-'S1S2;S3S4 — -L siS2;S3S4 ~ ^SiS2;S3Si, yKjOV ) 

J - f(i){o,-i} _ p (CAQ) 

^S1S2;S3S4 ~ SiS2;S3S4 -'S1S2;S3S47 I^V^IVJ^ 

L,,,,;.3,, = rilH2,0}^ ^ {U/16){S^S^ + 52S4 + 2Sl52S3S4)<5^ (C41) 

M - r(i)<^° "2'' — L (CA2) 

J"siS2;s3S4 ~ siS2;s3S4 -^siS2;s3S4) V^^^y 



p = r(i){3,-i} ^ _34 CC43') 

^ S1S2;S3S4 SiS2;S3S4 ^^S1S2;S3S41 V^^^^ 

n - p(i){i,-3} _ p CC441 

^;SlS2;S3S4 ~ -*- SiS2;S35-4 ~ -^ SlS2 ;S3S4 J V^^^^ 

i?,,,,;,3,, = r(l»t^;<?4^ - -(C//8)S1S2S354,52, (C45) 

'JsiS2;S3S4 ~ -■- SlS2;S3S4 ~ -'^SiS2;S3S4 (^v^^iu; 

The second-order amplitude is derived from Eq. (j4.22p with n — 2 

-'- SlS2;S4S3(,k, — k; p, — p) = —L y ^ \U s-4S2:,SsU SS;S3S4 + -tisiS2;SsJss;S3S4 

s 

-i-W V A- F n -\- J H ip*«Pp~«Sk I [w w 1/7 p ip2iep -2jek 

^^ L'^siS2;ss^ss;S3S4 ^^-'-siS2;ss^-'ss;s3S4 ^^<-'siS2:ss-'-'ss:S3S4j^ *-- ^^ L^*^siS2;ss'^^ss;S3S4 ^^^-'siS2;ss-*ss;s3S4j'^ ^ 

-1-4 T/ pidppi-^k IT/ D -iflp -iffk I FT/ p I p W 1 2jep -iSk 

^^ ^^siS2 ;ss*^ss;S3S4'^ ^ i^^siS2;ss-'--'ss;S3S4*^ "^ ^^L'^siS2;ss-'-ss;S3S4^^-'-siS2;ss'*ss;s3S4j'^ "^ 



„2i0„ 



S1S2 ;ss *^ ss;S3S4 •-■ ^ 1 ^ siS2;ss-'--'ss;S3S4^ ^ r [ f sis2;ss-^ SS;S3S4 I -^ siS2;ss ' * SS;S3S4j 

[^siS2;ss^ss;S3S4 I *-^siS2;ss ^ss;S3S4j^ ^ ^" [-^^5152 ;ss ^ss;S3S4 ~r f^siS2;ss-'^ss;S3S4j 

[ ^SiS2;SS<^SS;S3S4 T^ ^ S\S2\SS^ SS\S3S4\^ ^" [-'^SiS2;SS-^SS;S3S4 ~r tysj^s2;SS-'^SS;S3S4 ~r J-/ S\S2\SS*^* SS\S3S4\^ 

■\W M +n J -I- Af U ]p~'^^^^~i-V P p3*«Pp-«»k 

l'^*' SiS2',SS-'-^-^SS;S3S4 ^^ ^-^ SiS2',SS^ SS;S3S4 ^^ ■^*-^SiS2:SS^ ss;ssS4\^ T^ * siS2',ss-^ ss;ssS4^ ^ 

O V" p»^Pp~3iek I Ty n p4»ep -2i6lk I c w p^iSp -4i0kl 1 /n('x3\ ('C47'l 



- — Ws,s,-s,sAL) = Ws,sr,ss{L)Wss;s,sAL) + G s,srML)Fss;s,sAL) . (C50) 
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where L = (m/27r) ln(A/T). The second-order amphtude contains the same combinations of the harmonics e^^^p 
and e™^^ as the first-order amphtude, which proves the group property. The RG flow equations are obtained by 
replacing the left-hand side of Eq. (JC47P by the bare amplitude, letting the coefficients f/ ... 5 to depend on L, and 
differentiating with respect to L: 

rr^siS2:SiS4(^) = UsiS2:ss{L)Uss:s^Si{L) + Hs-^s2:ss{L)Jss-s3S4{L), (C48) 

dL 

~-TY'^siS2;S3Si(L) = Vs-^s2\ss{L)Vss-S3Si{L) + Fs-^s2;ss{L)G ss-S3Si{L) + J siS2;ss{L)H ss-s:-iSi{F) ■, (C49) 

d_ 

dZ' 

-^As,s2;S3sAL) = As,s2-ML)Vss-S3sAL), (C51) 

--Tj-BsiS2;s3Si{L) = Vs-^s2;ssiL)Bss;s3Si{L), (C52) 

— --JYFsiS2;S3S4{L) = VsiS2;Ss{L)Fss;S3S4{L) + i^si S2 ;SS (-^)W^SS;S3 S4 (^) J (C53) 

—^Gs,s2;S3sAL) = Ws,s2-ML)Gss;S3sAL)+Gs,s2-ML)Vss;S3sAL), (C54) 

^~rr-^'>-l-'^2-S3S4{F) = Hs^S2;Ss{L)Vss-S3S4\L) + f/^ j 53 ;ss (i)-ffsS;S3 S4 (-t') , (C55) 

dL 

~'Tj'JsiS2;S3S4{L) = Vs-i^s2;ss[L)Jss;S3S4(Lj + JsiS2;ss[L)Uss-.S3S4[L), (Cod) 

^'Tj'FsiS2\S3S4{L) = Hs-^s2;ss{L)Fss-s3S4{L) + [/g j gj ;ss (I/)iss;S3S4 (-^) + Ls-i^s2;ss{L)Wss;S3S4{L), (C57) 



Ms,S2;S3S4iL)^Ws,S2;Ss{L)Mss;S3S4{L)+Gs,S2-ML)Jss;S3S4{L)+Ms,S2;Ss{L)Uss;S3S4{L), (C58) 



_d_ 
dL 

~'Tj'FsiS2;S3S4{L) = VsiS2;Ss{L)Pss;S3S4{L), (C59) 

-^Qs,S2;S3S4{L) = Qs,S2-ML)Vs,;S3S4{L), (C6O) 

A 

dl 



— -J^RsiS2;S3S4iF) — WsiS2lSsiL)RssiS3S4iL), (C61) 



dL 



SsiS2;S3S4{F) — Ss-iS2\Ss{L)Wss-S3S4{Lj, (C62) 



where summation over the repeated index s is implied. The initial conditions are given by X^^ ,52;s3,s4(0) — Xs-^^S2;S3,S4 
with X — U . . . S. Since it is very difficult to solve this system of differential equations analytically, a new approach 
is required. In what follows, we will determine U{L), . . . , S{L) for a few lowest orders in the "RG time" L and then 
make a guess for a form of an arbitrary-order term. The RG equations will then provide a necessary check as to 
whether our guess, based on the perturbative calculation, gives a correct answer. 

A few lowest order amplitudes can be derived perturbatively from Eq. (j4.22p . copied here for the reader's convenience 

ril;a4S3(k,-k;p,-p)(L) = -i5]y'"''^ri^;,lL(k,-k;i,-i)rW,^,3(i,-i;p,-p) (C63) 

with j > 2 standing for order of the perturbation theory. Since the scattering amplitudes depend on angles 9^ and 
6p, they can be decomposed order by order into the Fourier series 

oo 

ri^^),^,3,Jk,~k;p,-p)= Y^ r(^ll".3"le™'''e"''N (C64) 

rn,n— — oo 

where the coefficients in front of e^"^^pe^"^^ are determined using the orthogonahty property 

m^Ts^t (k, -k; P, ^P) = l'^ ^ fj ^m^;s4S3 (k, -k; p, -p)e-™^pe-»^^ (C65) 
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Resumming the coefficients of g*™Spe»"^k ^q infinite order (with m and n being fixed) 

oo 

ririlr:^:^ (k, -k; p, -p) = Y. ^ii&t (k, -k; P, -p) (Cee) 

we can find the renormalized amphtudes. For each combination of partial harmonics, which occurs to second order 
in the magnetic field, we derive explicitly the scattering amplitudes up to seventh order in the Cooper channel 
renormalization parameter UL and then make a guess for general j-th order amplitude. The final result is obtained 
by resumming these amplitudes to infinite order and then substituted into the RG flow equations to check the 
correctness of our guess. In all cases, the guess turns out to be correct. All nonzero RG charges as well as their large 
L limits are listed below. We begin with the n = m = and n = m = 1 harmonics, given by 

Uss-ML) - JYUlL + J^' ^ ¥ ^' + ^(^""' ^^' ^^^^^ 

Us-s.^-UL) 4^ - f ^^^ - Oil.-^ n (C69) 

t/,_,;,_,(L) = ^-^ + ^-^6^ = Oiln-' T), (C70) 



2 2 + UL 8 1 + UL 
_ U 1 



f/.(±T;TT) W = ^TT^ = ^(1^"' T), (C71) 



V,s-ss{L) = ^- ^——7^6^ = 0(ln-^ T), (C72) 

Vs,-,-s(L) = i + ^——^S^ = 0(ln-^T), (C73) 

ss, s s\ J 2 1 + UL 4 {1 + ULf V ;: \ > 

F._,;_,.(L) = -^ + ^(1 + UL)5^ = -^ + 0(C/2), (C74) 

Vs-s;s-s{L) = I - ^(1 + UL)5^ = ^ + 0{U'). (C75) 

An important remark should be made at this point: in addition to amplitudes which flow either to zero or to finite 
values at low temperatures, there are also amplitudes which grow logarithmically at low temperatures, namely, the 
amplitudes in Eqs. (|C74|) and (jC75p . This peculiar feature, which occurs only in the presence of both the SOI and 
in-plane magnetic field, may indicate a phase transition below certain field-dependent temperature or it may be an 
artifact of the expansion to lowest order in S^. In the derivation of the spin susceptibility that follows in App. IC 2bl 
we assume that the electron gas is far above the temperature below which the instability becomes important, i.e., 
that UL <^ 1/(5^, so that the effect of the instability can be neglected but the nonperturbative regime of Cooper 
renormalization, where 1 <C UL <C l/S"^, can still be accessed. 
The remaining harmonics are 

W,sML)^Uss-ssiL), (C76) 

W,s--s-s{L) = Uss--.-s{L), (C77) 

W,.s;^ss{L) - Us-s--ss{L), (C78) 

W,.s;s-s{L) = Us-s-s-s{L), (C79) 

W^.(±t;tt)(^) = -Ua(±T;TT)iL), (C80) 

A^s-ML) = -^,;_,_,(i) = ^-^-^-^S^ = 0(ln-i T), (C81) 

A,_s;-ss{L) = -A,_,;,_,(i) = -^S^ (C82) 

lb 

BsiS2;S3Si[L) — As-^s2-s3Si[Lj. (C83) 
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F±^.,^^{L) = ^^i6, (C84) 

i^T±;TT(i) = J^S, (C85) 

F^^.,i^{L) = -^^-^z6 = 0(ln-'T), (C86) 

G±^;ttW = JY^^*<5 = 0(ln"'r), (CSS) 

G^T;±T(i) - J*'^, (C90) 

G^^;^±(L) = -^zJ, (C91) 

HsiS2;S3S4,{L) — Gs-iS2;SsS4{L) (C92) 

JsiS2lS3S4{L) = FsiS2:S3S4{L) (C93) 

^--(^) = ? (2tW + ^'^ ^ ^'^ + ^^^^" '^^' ^^''^ 

i.;-.-. W - ^ (^^ - f ^^ - -J^' + ^(In- T), (C95) 

' ^ ' 8 (l + t/L)(2 + t/i)2 V ^ V ; 

^^---^(^) = ¥(m7L)(2TT7LF ^ ^ n (C97) 

Mss;sAL)=Lss;ss{L), (ClOO) 

M,,;_,_,(L) = L,,;_,_,(L), (ClOl) 

Ms-s;-ss{L) = L,_,;_,,(L), (C102) 

M,^s;s-s{L) = L,_,;,_,(L), (C103) 

M±^;^^(i) - L^±;tt(L) = -^^^-1^52 = 0(ln-i T), (C104) 

Mtt;±^(L) = L^^,^±{L) = -^A±RL_5^ = o(in-i T), (C105) 

"siS2;S3S4(i) = QsiS2;S3S4 l-'^) — ^'■^'^SiS2:S3S4,\L), (C106) 

Rss,ss{L) - Rss,-s-s{L) = Rs-s;-ss{L) = Rs-s;s-s{L) = -JJTUL^^ " ^^^''"' ^^' ^^^°^^ 

i?.(±T;TT)W = J2TC7I'^' = 0(ln-'T), (C108) 
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SsiS2;S3Si{L) = RsiS2;S3SiiL). (C109) 

It can be readily verified that all the amplitudes satisfy RG equations (|C48p - (jC62[) with initial conditions (IC32l) - (jC46l 
up to 0{S^) accuracy. 

Finally, the renormalized scattering amplitude is given by 



lSiS2;S3S4 



A,,,,;,3,,(L)e'''''e*'''^ + B,,,,;,3,,(L)e"'''''e^*'''^ + ^.,,,;,3,,(L)e2'''>'e-'''>' + G,,,,;.3.,(L)e*''-e- 



(Clio) 

b. Renormalization of Xxx 

As for the transverse-field case, the free energy for the in-plane magnetic field is found by replacing the bare 
interaction U in Eq. p. 381) by the renormalized vertex T 

s^..^-\ r^TY^T. /°"§r,,,,,,3,,(k,-k;-k,k)r.3.,,.,,,(~k,k;k,-k)n+':-n-'^,-, (cm) 

where 11^^,, ^ given by Eq. p.39p depends on the direction of the electron momentum with respect to the magnetic 
field. 

A general formula for S'^xx is very complicated; however, in the regime of strong Cooper renormalization, i.e. for 
1 ^ UL ^ 1/(5^, there are only a few partial amplitudes which survive the downward renormalization. Keeping only 
these partial amplitudes in F, we obtain for the thermodynamic potential 

6Exx = -^fJ ^tY^ J '-^mll^nzl^ + n+1-n;>l- - 2ng) + ml] 



A' U 



rE/?[no(n-+ + n+_-2no) + n2], (cii2) 



a^kl 8 ^ 

where the angular dependence of the polarization bubbles in the second line was neglected because of an overall factor 
A^ originating from the scattering amplitudes, and the angular integral in those terms was readily performed. On 
the other hand, the field dependence in the first term is exclusively due to the bubbles, hence the angular integration 
has to be carried out last. The integrals over the momentum and frequency yield 



^^xx 



T^ fmUV\ /-^"rfek^^/ Ak, +A_k, \ ^ ^ A^ ^^f\a\k 



87rv|, V 47r 



27r V T ) a^Pp \ T 



(C113) 



Expanding J-{x) ss a;'^/3 for x ^ 1 and differentiating with respect to the field twice, we obtain for the nonanalytic 
part of the spin susceptibility 

1 (mUV\a\kp 
5xxx^-^X.[^-^) -^. (C114) 

Somewhat unexpectedly, the fully renormalized result (jC114[) coincides with the leading (first) term in the second- 
order result (|3.51l) . The formally subleading but T-dependent T/2Ef term in Eq. (|3.5ip does not show up in the fully 
renormalized result, which implies that, at best, it is of order T JUL ex T/ InT for large but finite UL. Hence follows 
the result for &Xxx presented in the main text, Eq. (|4.39p . 
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